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ABSTRACT 


A preliminary,  mathematical  paper  [4]  set  forth  the  theoretical 
basis  of  permanent  memory  behavior  in  a restricted  class  of  polymeric 
solids.  The  present  report  expands  considerably  upon  [4]  and  sets 
forth  a 3-dimensional  theory  of  permanent  memory  effects.  Metatropic 
Materials,  applicable  to  heterogeneous  solids. 

The  change  in  the  value  of  the  various  moduli  as  well  as  the 
changes  in  the  isotropy  group,  i.e.  crystal  class,  of  the  stress 
response  caused  by  a past  sequence  of  strain  histories  is  given  a 
general  structure. 

The  nature  and  sign  of  the  residual  stresses,  permanent  sets,  and 
the  isotropy  group  are  deduced  from  an  operational  definition  of 
metatropic  effects.  These  include  damage  effects,  strain  hardening, 
and  aging. 
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CHAPTER  I : OVERVIEW 


1.1  Introduction 

Various  theories  of  damage,  permanent  set,  strain  hardening  and 
residual  stress  have  appeared  in  the  literature  over  the  past  100  years. 
Most  notable  of  the  earlier  works  were  those  of  Lord  Kelvin  [1]  and 
E.  Warburg  [2].  Kelvin  made  the  observation  that  the  present  state  of 
a material  element  depends  at  any  instant  on  its  previous  states  as 
well  as  on  the  stimuli  which  obtain  at  the  instant.  He  was,  in  this 
context,  referring  to  permanent,  irreversible  changes  in  the  material 
and  not  just  to  the  viscoelastic  or  fading  memory  aspects  so  much  in 
vogue  today. 

Warburg,  in  a series  of  experiments,  demonstrated  two  very  impor- 
tant phenomena  particularly  applicable  to  observations  on  filled  and 
fibre  reinforced  polymers; 

(1)  A material  after  moderate  shear  straining  exhibits  a decrease 
in  its  tensile  modulus  and  an  increase  in  the  value  of  its 
compressive  modulus  and 

(2)  An  initially  isotropic  material  after  shear  straining 
becomes  essentially  orthorhombic  in  its  subsequent  crystal 
class  mechanical  response. 

Warburg  was  examining  three  dimensional  permanent  memory  effects, 
and  observing  changes  in  both  the  values  of  the  various  moduli  as  well 
as  changes  in  the  isotropy  group  of  the  material  elements. 


An  excellent  review  by  Bartenev  and  Zuyez  [3]  points  out  rather 

strongly  that  " creep  and  incipient  failure  proceed  simultaneously- 

The  elementary  act  of  failure  precedes  the  elementary  act  of  creep". 
They  are  referring  to  what  we  herein  term  Permanent  Memory  effects. 
Bartenev  also  stresses  the  great  lack  of  true  three-dimensional  theory 
applicable  to  polymers. 

It  is  the  purpose  of  this  paper  to  determine  from  a phenomenolo- 
gical viewpoint  the  assumptions  necessary  to  formulate  a three-dimen- 
sional framework  for  Permanent  Memory  effects  and  to  examine  and  set 
forth  the  necessary  mechanical  consequences  of  such  a generalized 
theory  with  respect  to  moduli,  permanent  set,  residual  stresses,  and 
changes  in  the  isotropy.  group  of  the  material  element. 

1.2  Metatropic  Materials 

The  senior  author  [4]  set  forth  an  initial  attempt  with  the  above 
mentioned  purpose  in  1974.  At  that  time,  the  effect  of  prior  strain 
histories  on  the  istoropy  group,  i.e.,  crystal  class,  of  the  filled 
polymers  tested  was  not  included  nor  were  the  conclusions  relative  to 
permanent  set  and  residual  stresses  fully  formulated.  Damage  was  the 
term  then  used  to  label  these  permanent  memory  phenomena  and  included 
the  rather  unsatisfactory  concept  of  positive  and  negative  damage  ef- 
fects. Other  authors  [5]  attempting  to  correct  the  nomenclature 
problem  coined  the  equally  unsatisfactory  term  "divagation"  meaning 
material  whose  properties  wander  about. 
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Herein  we  propose  the  term  Metatropic  Materials  defined  as  "mater- 


ials which  change  or  tend  to  change  their  properties  in  a specified 
manner  (meta-)  and  which  produce  or  tend  to  produce  a more  highly 
organized  state  (-tropic)  in  response  to  the  sequence  of  past  stimuli". 

While  metamorphic  might  seem  a similarly  suitable  term,  its  his- 
toric use  in  geology  could  confuse  the  issue. 

The  absence  of  the  metatropic  effects  results  in  Equitropic 
materials  whose  isotropy  group  is  invariant  after  any  imposed  strain 
history.  This  class  of  materials  includes  the  usual  linear  or  non- 
linear viscoelastic  solids  and  fluids  as  well  as  Noll’s  [6]  semi- 
elastic materials. 


1.3  The  Metatropic  Function 

We  denote  the  metatropic  function,  7YI,  and  normalize  its  value  to 
unity  in  the  virgin  state  of  the  material. 

It  can  be  readily  shown  that  the  metatropic  function,  whether 
signifying  a damaged  material,  7/1 < 1 or  a strain-hardened  material, 
m>i,  can  only  be  determined  in  general  if 

(1)  The  two  samples  being  compared  are  in  the  same  configuration 
and 

(2)  Both  are  completely  relaxed  with  respect  to  creep  or  stress 
relaxation  effects. 

As  an  example,  consider  two  samples  of  an  ideal  linearly  visco- 
elastic material  wherein  777,  =1  for  all  states.  Let  <t,t  denote  the 
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state  of  each  sample  respectively  and  let  them  each  be  in  identical 
configurations,  G. 

Noll  [6]  defines  equality  of  state  operationally,  that  is, 
given  two  material  elements  in  the  same  configuration,  G,  subject  to 
a deformation  (strain)  process,  P,  if  the  observed  tensorial  stress 
outputs,  S,  are  such  that 

S(a,P)  = S(t,P) 

for  every  process,  P,  then  we  say  that  the  states  are  equal,  i.e., 
o = t.  If  inequality  holds  for  any  P,  then  the  states  o,t  are  unequal. 

If  the  material  element  in  the  state  a is  one  wherein  the  effects 
of  a previous  strain  history  have  not  fully  relaxed  and,  say,  t 
represents  the  virgin,  relaxed  state  then  obviously  the  stress  responses 
will  not  be  equal  even  for  our  linearly  viscoelastic  model. 

Thus,  in  order  to  determine  whether  the  states  a and  t differ 
with  a permanent  metatropic  change,  all  comparisons  must  be  carried 
out  from  fully  relaxed  states  in  identical  configurations. 

We  wish  to  associate  the  primitive  concept  of  damage  with  some 
sort  of  weakening  effect  and  strain-hardening  with  some  sort  of  strength- 
ening effect.  The  inequality  of  Fitzgerald  [4]  embodies  these  physi- 
cally desirable  characteristics  as  follows: 

Definition  1 

For  states  a, r of  two  fully  relaxed  material  elements  in  identical 
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configurations,  G,  we  say  that  the  metatropic  function  171/  (o)_<  7)%/ ( t) 
if  and  only  if  the  stress  response  functional  S(a,P)£  S(t,P)  for  every 
deformation  process  which  starts  from  the  configuration,  G. 

Mathematically,  we  have  defined  a partial  ordering  on  Ht  through 
an  isomorphism.  States  which  obey  D1  are  called  comparable. 

States  which  do  not  obey  D1  are  called  incomparable  and  denoted  by 


For  uniaxial  tests,  all  states  are  comparable.  For  general  three- 
dimensional  strain  histories,  however,  the  tensorial  nature  of  stress 
can  produce  incomparable  states. 

1.  A Tensorial  Character  of  m, 

Let  us  examine  the  tensorial  character  of  m from  an  heuristic 
point  of  view.  For  fully  relaxed  states  in  identical  configurations 
subject  to  slow,  infinitesimal  strains,  e rather  than  general  de- 
formation processes,  P,  the  linear  elastic  relation  is 

ij  ijkl  kl 

wherein  is  a fourth  order  tensor  in  general, 

ijkl 

The  metatropic  function,  m , induces  a change  in  the  material 
coefficient  values  of  the  elasticities,  C as  well  as  possibly  changing 
the  isotropy  group  or  crystal  class  of  C.  Thus  the  linear  transformation 
from  an  initial  elasticity  to  a changed  elasticity  is  ex- 

pressed by 
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•->-  • *• 


•J 


C11kl  = ^ c(1) 

ijkl  ijklmnop  mnop 

with  in  an  eighth  order  tensor. 

Taking  into  account  the  usual  symmetries  of  C,  this  linear  approx- 
imation to  TYl  yields  an  expression  requiring,  in  general,  the  eva- 
luation of  231  constants.  The  impracticability  of  experimentally 
determining  these  constants  is  obvious. 

Since  experimental  results  on  filled  amorphous  polymers  show  a 
certain  "smearing"  with  respect  to  crystal  class  axes  after  a damaging 
history,  it  becomes  practicable  to  restrict  the  form  of  m to  no  more 
than  that  of  a second  order,  symmetric  tensor  as  was  done  in  [A].  This 
restriction,  however,  is  only  rigorously  applicable  to  strain  histories 
whose  principal  directions  do  not  vary  with  time  even  though  their 
principal  values  vary  arbitrarily. 


The  restriction  then  implies  that,  for  example,  strain  history 
damage  effects  can  change  an  initially  isotropic  material  into  a 
transversely  isotropic  one  and  to,  at  most,  an  orthorhombic  material. 

1.5  Physical  Consequences  of  "WL 


The  use  of  the  inequality  in  D1  and  a few  added  mathematical 
assertions  lends  to  the  following  set  of  results  : 

(1)  Prior  straining  causes  a change  in  the  stress  response 
moduli 


* 


(2)  Damage  produces  compressive  stresses  in  the  unity  configura- 
tion and  an  "elongation"  type  of  permanent  set  in  the  stress- 
free  configuration 

(3)  "Strain-hardening"  and  aging  produce  tensile  stresses  in  the 
unity  configuration  and  a "compressive"  type  of  permanent 
set  in  the  stress-free  configuration 

(4)  Prior  straining  changes  an  initially  isotropic  material 
into  an  aelotropic  material 

(5)  Changes  in  the  structural  component  of  a material  cause 
"elastomorphic"  behavior  changes 

(6)  Changes  in  the  viscous  component  of  a material  cause 
"viscomorphic"  behavior 

(7)  A viscoelastic  range  can  exist  for  most  filled  materials. 


1.6  Summary  of  Work  to  Date 

In  order  to  provide  an  overview  of  the  development  of  the  ideas 
and  methods  herein  presented,  the  next  several  pages  constitute  a 
summary  of  the  effort  with  pertinent  references. 

Damage  is  defined  as  a modulus  reduction  caused  by  any  micro- 
structural  weakening  in  a material  arising  from  prior  deformation 
histories  (damage  may  be  inherent  in  the  material  as  a result  of 
processing,  then  the  "virgin"  state  may  be  inferred  by  a limit  test- 
ing process  but  the  "virgin"  state  is  not  accessible) . 
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Consider  two  material  samples  to  be  compared  for  (tensorial) 
damage  with  S,  the  stress  tensor;  n an  outward  acting  unit  normal; 

D,  the  damage  tensor;  P,  a deformation  process,  and  G,  the  confi- 
guration in  a relaxed  state,  we  then  say  [7]  the  degree  of  damage 
D2  < Di  if  and  only  if  for  every  deformation  process  P starting  from 
identical  relaxed  configurations  G,  the  magnitudes, of  the  projec- 
tion of  the  surface  tractions,  S-,  projected  onto  the  exterior  unit 
normal  n,  are  related  by  a dual  isomorphism* 


(Si  n)  *n  <_  (S2  n)  .n 


(1-1) 


with  strict  inequality  holding  in  (1-1)  for  at  least  one  P and  one 
n.  Equality  in  (1-1)  for  every  P,  n means  equality  of  damages 
D1  = 1*2 • (Note:  D2  < Di  - (Sin)*n£  (S^i)  *n,  i.e. , leads  to  a 

"weaker"  material,  lesser  modulus.) 

* Note:  the  damage  tensor  is  generally  of  second  order  for  practical 

reasons  and  is  in  an  inverse  relation  to  T[Q  • 


1.6.1  Comparability  of  Damage 

Since  D is  a symmetric  tensor,  not  all  D's  may  be  comparable 
under  the  conditions  of  Equation  (l-l). 


We  say  that  damage  states  obeying  (l-l)are  comparable.  Other- 
wise, they  are  incomparable. 
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(Example:  one-directlonal  damage  may  have  a weakening  by  (1-1)  in 

that  direction  compared,  to  say,  spherical  (isotropic)  damage  in 
another  specimen.  But  the  spherical  damage  specimen  will  be  weaker 
in  a transverse  direction,  hence,  the  conditions  of  (l-l)are  violated 
and  the  two  damage  states  are  incomparable,  i.e.,  effectively  consist 
of  two  different  materials) . 

Thus  we  achieve  a partial  ordering  on  damage,  D,  so  that  our 
topology  is  a space  of  chains. 

This  implies  (Zorn's  Lemma)  that  a maximal  element  exists,  i.e., 
failure  when  from  (1-1) 

(Sn)  *n  0 

(Either  no  stress  can  deform  specimen,  broken,  or  only  compressive 
stresses  can  be  sustained,  i.e.,  "sticking"  it  back  together.) 

We  postulate  a virgin  state  which  is  the  minimal  element.  Since 
(l-l)is  a scalar  number  ordered,  on  a given  chain,  by  damage  D,  we 
thus  have  a material  ordering  on  D on  a chain  which  can  be  normalized, 
say  from  0 to  1 for  comparable  damages.  [7] 


1.6.2  Semi-elastic  Range  (A  Visco-elastic  Range) 

So  long  as  a material  is  subject  to  deformation  processes  in 
its  semi-elastic  (visco-elastic)  range,  no  further  damage  ensues 
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and  Its  relaxation  modulus  Remains  unchanged. 


The  term  semi-elastic  range  Is  from  [ 7 ] and  visco-elastic  range 
is  defined  and  demonstrated  in  [ 8 ] and  [ 9]. 

Note:  For  some  materials  a semi-elastic  range  may  not  exist,  i.e., 

any  deformation  or  motion  may  lead  to  increased  damage. 


1.6.3  Residual  Stress 

From  [ 4]  it  is  shown  that  if  a damaged  material  is  forced  back 
into  its  unity  (i.e.,  virgin)  configuration,  then  residual  compressive 
stresses  will  result.  Proof:  from  (1-1)  Sjn-n  <_  S2n*n  if  D2  < Dj. 

Consider  material  2 as  virgin.  Then  with  both  samples  in  original 
undistorted  virgin  configuration,  the  stress  needed  to  maintain  this 
configuration  on  the  virgin  sample  is  clearly  zero,  so 

Sin*n  <_  0. 

That  is,  any  residual  stresses  on  the  damaged  sample  are  zero  or 
compressive. 


1.6.4  Residual  Stress  Compressive  Only 


In  [io  ] it 
residual  stress 


is  then  proven  that  they  cannot  be  zero,  hence,  the 
field  on  a damaged  sample  to  force  it  into  a (relaxed) 
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"virgin"  configuration  are  somewhere  compressive  and  possibly  some- 
where zero  but  are  nowhere  tensile. 

1.6. 5 Permanent  Set 

Postulating  certain  injectivity  properties  for  the  relaxed  con- 
figuration defines  various  classes  of  materials  [ 7,10].  These 
classes  are  shown  in  Figure  1,  attached. 

It  is  then  proven  in  [10]  that  a damaged  material  must  have  per- 
manent set  in  the  relaxed,  stress-free  configuration.  It  is  further 
proven  that  the  permanent  set  configuration  (3-D)  is  ordered  by 
the  partial  ordering  on  damage,  that  is,  given  a chain  of  comparable 
damage  states,  the  permanent  set  will  be  uniquely  determined. 

The  inverse  is,  given  a series  of  permanent  sets  caused  by 
comparable  damage  states,  the  set  uniquely  determines  the  degree  of 
damage. 

1.6.6  Permanent  Set  "Extensile"  in  Nature 

For  a relaxed  configuration  it  is  shown  in  [12]  that  the  stress 
is  related  to  the  deformation  (strain  measure)  by  a pseudo-linear 
relation 


S = GC 


(1-2) 


with  S,  stress;  G,  a damage  and  "strain"  dependent  modulus  (4th  order 
tensor)  and  C,  the  strain  measure.  Equation  (2-2) is  the  1st  order 


■ 
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The  presence  of  the  fourth  integral  leads  to  an  apparent 
"stiffening"  of  the  modulus  under  active  damaging,  i.e.,  the  material 
stress  response  shows  a steepening  with  respect  to  strain  when  re- 
approaching the  "virgin"  curve  as  in  the  Mullins'  effect  repeated 
loadings . 

A computerized  characterization  method  using  (1-3)  has  been 
developed  and  is  being  used  by  United  Technology  Corporation  on  an 
investigative  basis.  The  general  approach  is  in  [12,13]. 

Other  characterizations  such  as  those  of  Farris  and  Schapery  are 
sub-classes  of  the  above  Equation  (1-3)  as  shown  in  [14]. 

1.6.8  Thermo-mechanical  Coupling 

Lee  [15]  has  derived  an  equation 

S = GjC  + /G2(D,0,C,t-x)  C(x)dx 
+ /G3(D,e,C,t-T)  D(x)dx 

+ /GH(D,9,C,t-x)  0(x)dx  + G5D  (1-4) 

where  6 is  the  temperature  excursion  from  a rest  state.  For  the 
case  of  D = 0,  he  has  found  numerical  values  for  G4  and  has  success- 
fully accounted  for  simultaneous  cooling  and  straining  results  in 
damaged  materials. 

1.6.9  Chemical  Aging 

Hufferd,  using  a form  of  (1-4)  but  with  H substituting  for  0 
is  using  it  with  Thiokol  (John  Bennett)  to  account  for  chemical 
aging  effects. 
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1.6.10  Hierarchy 

The  hierarchy  of  materials  is  given  in  Figure  1 . 

Semi-elastic  means  general  linear  or  nonlinear  visco-elastic. 
Elastic  means  linear  or  nonlinear. 

Mullins  elastic  covers  elastic  materials  with  Mullins  effects 
and/or  hysteresis. 

All  are  sub-classes  of  my  PMSE  type,  permanent  memory  with  semi- 
elastic range. 

(I  adopt  the  word  semi-elastic  from  Walter  Noll,  whose  methods 
I have  used  to  frame  my  PMSE  concepts.) 

This  figure  will  appear  in  [10  ] . 


1.6.11 


Metatropic  Materials 


Observing  that  prior  damage  producing  straining  processes  effect 
not  only  a change  in  the  value  of  the  response  functional  but  also  pro- 
duce a change  in  the  isotropy  group  of  the  material,  the  term 
metatropic  was  coined  by  the  author  to  describe  these  combined  effects. 

Metatropic  Materials  are  defined  as  "materials  which  change  or 
tend  to  change  their  properties  in  a specified  manner  (meta-)  and  which 
produce  or  tend  to  produce  a more  highly  organized  state  (-tropic)  in 
response  to  the  sequence  of  past  stimuli". 

The  above  is  presented  in  [16]. 
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1.5.12  Finite  Deformation  Isotropy  Groups 


A parallel  development  has  been  the  observation  (rediscovery?) 
that  materials  may  have  higher  order  symmetries  which  appear  isotropic 
under  small  strain.  For  example,  a material  which  has  transverse  octa- 
gonal symmetry  will  exhibit  identical  response  at  0°,  45°,  90°, 
at  infinitesimal  strain  and  thus  will  be  isotropic. 

However,  at  large  strains,  the  octagonal  symmetry  will  be  ob- 
served. In  general,  for  each  higher  polynomial  in  strain  that  is  con- 
tained in  the  strain  energy  function  for  a material,  two  additional 
symmetry  groups  can  be  observed. 

These  results  have  application  to  fibre  reinforced  materials 
undergoing  large  strains.  This  work  is  covered  in  [ 17  ]. 

A loosely  related  paper  evolving  from  the  ideas  developed  herein 
is  [ 18]. 
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NOMENCLATURE,  CHAPTER  II 


„ JO)  (1) 

c , c * c 


Modulus  tensors 
Damage  tensor 

Eulerian  damage  measure  (2.8) 

Damage  in  a linear  theory  (2.10) 

Green  strain 

Green  strain  referred  to  damaged  state 
Eulerian  strain 
Infinitesimal  strain 

Deformation  gradient 
Inverse  deformation  gradient 

Metric  tensor  in  undeformed  state 
Metric  tensor  in  damaged  state 
Metric  tensor  in  deformed  state 

Load  history 

Stress  in  a linear  theory 

Symmetric  Kirchhoff  stress  tensor 
Symmetric  Kirchhoff  stress  in  damaged  state 
Cauchy  stress  tensor 

Cartesian  coordinates  in  undeformed  state 
Cartesian  coordinates  in  deformed  state 


CHAPTER  II:  RECENT  WORK 


2.1  This  report  presents  recent  work  on  permanent  memory  materials. 
Permanent  memory  phenomena  were  first  observed  in  uniaxial  tests  of  highly 
filled  solid  propellants.  When  specimens  were  strained,  brought  to  rest, 
and  then  retested,  a reduction  in  moduli  was  observed.  Moreover,  an 
increase  in  the  applied  strain  led  to  a larger  modulus  reduction.  The 
term  "permanent  memory"  was  introduced  to  distinguish  this  behavior  from 
the  "fading  memory"  behavior  of  viscoelastic  materials.  The  modulus 
reduction  caused  by  strain  is  essentially  irrecoverable;  it  is 
definitely  not  a transient  phenomenon. 


In  the  case  of  propellants,  straining  causes  a degradation  of 
material  properties.  Therefore,  the  change  was  called  "damage".  Sub- 
sequently, permanent  memory  properties  were  found  to  be  fairly  common ; 
they  do  not  always  lead  to  degradation,  but  may  result  in  improvement. 
Thus,  by  permanent  memory  we  now  mean  any  permanent  change  in  properties 
due  to  strain  history.  The  more  specific  term  "damage"  is  retained  here; 
nevertheless,  it  will  be  clear  that  our  theory,  suitably  modified,  can 
also  model  improvements  in  material  properties. 


To  formulate  the  problem  discussed  here,  suppose  that  a material 
is  in  some  initial  state  in  which  its  mechanical  properties  are  known. 
Next,  let  it  be  subjected  to  arbitrary  loads  and  then  brought  to  rest; 
our  problem  is  to  predict  the  subsequent  mechanical  behavior.  Because 
of  permanent  memory  effects,  it  will  not  be  the  same  as  it  was  before 
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Our  treatment  of  this  problem  is  based  on  certain  key  ideas  from 
previous  work.  The  first  is  the  idea  of  a viscoelastic  range.  This 
term  denotes  strain  histories  to  which  a permanent  memory  material 
responds  viscoelastically.  Thus,  in  a viscoelastic  range,  no  further 
change  in  material  properties  occurs.  This  idea  simplifies  our  problem 
considerably.  It  now  becomes  a problem  of  relating  viscoelastic 
characterizations  before  and  after  damage.  Second,  we  use  the  suggestion 
of  Fitzgerald  [1]  that  damage  be  modeled  by  a second  rank  tensor.  Our 
work  actually  sheds  some  light  on  when  this  model  is  valid.  It  leads 
to  a definite  mathematical  problem  which  is  ultimately  solved  by  Rivlin's 
theory  of  tensor  valued  functions.  Third,  we  use  the  experimental 
observation  that  damage  is  caused  by  tensile  strains.  From  this  idea, 
we  can  derive  definite  expressions  relating  the  damage  tensor  to  the 
strain  history.  In  summary,  our  basic  premises  are: 

(1)  The  existence  of  a viscoelastic  range 

(2)  Damage  is  described  by  a second  tank  tensor. 

(3)  Damage  results  from  tensile  strains. 

In  section  2,  the  basic  problem  is  formulated  and  discussed  in 
greater  detail.  Section  4 introduces  the  tensorial  damage  measure. 

Damage  dependent  stress-strain  relations  are  developed  in  sections  4 and 
3.  Section  6 describes  the  qualitative  effects  which  this  theory  can 
model.  In  section  7,  the  stress-strain  relations  are  rederived  from 
Green's  point  of  view  by  using  a damage- dependent  strain  energy  function. 
In  section  8,  the  stress-strain  relations  are  written  in  Eulerian 
variables;  in  section  9 they  are  written  in  terms  of  the  damaged 
traction-free  state.  All  of  this  work  is  valid  for  finite  deformations. 
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In  section  10,  a linear  theory  for  permanent  memory  materials  is 
proposed. 

The  basic  results  are  equations  (5-2)  and  Q0-1) . Equation  (5-2) 
is  the  finite  damage  dependent  stress-strain  relation  for  an  initially 
isotropic  material.  The  assumption  and  simplifications  which  lead  to 
(5-2)  are  carefully  stated  in  sections  2 through  5.  Equation  (10-1)  is 
the  corresponding  formula  for  infinitesimal  strains.  According  to  this 
theory,  the  damaged  material  will  contain  permanent  strains  relative  to 
the  undamaged  material.  Also,  the  damaged  material  may  be  anisotropic. 
The  minimum  symmetry  possible  is  orthorhombic,  with  crystal  axes 
coincident  with  the  principal  axes  of  the  damage  tensor.  Section  6 
contains  a detailed  discussion  of  the  symmetry  properties  of  the  damaged 
material . 

Finally,  we  note  that  our  goal  has  been  to  derive  the  most  general 
equations  consistent  with  our  assumptions.  This  goal  is  in  keeping  with 
Rivlin's  approach  to  deriving  stress-strain  relations.  As  a result,  the 
theory  is  considerably  more  complex  than  necessary  for  materials  such  as 
propellants,  for  which  experimental  data  contains  a very  high  amount  of 
scatter.  However,  since  the  theory  is  general,  it  can  model  many 
different  qualitative  properties.  Once  one  determines  which  of  these 
properties  are  appropriate  for  a material,  characterization  only  requires 
choosing  a few  of  the  corresponding  terms  from  (5-2),  and  curve  fitting 
the  experimental  data. 

2.2  Suppose  that  a permanent  memory  material  has  the  stress-strain 
relation 
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This  formula  characterizes  the  material  in  the  virgin  state,  or  in  any 
fixed  reference  state.  The  material  is  loaded  arbitrarily  and  then  brought 
to  rest.  Because  of  permanent  memory  effects,  (2-1)  no  longer  describes 
the  material,  and  a new  characterization  is  needed: 
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In  (2-2),  the  term  reflects  the  possibility  of  permanent  set: 

nonzero  stresses  may  be  needed  to  restore  the  original  rest  configuration. 
The  problem  of  this  section  is  to  predict  formula  (2-2)  knowing  (2-1)  and 
the  intervening  loads.  Mathematically,  the  relationships 
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(2-3) 


are  sought,  where  L stands  for  the  damage  producting  loads. 


Posing  this  problem  is  justified  by  the  existence  of  a viscoelastic 
range  in  many  permanent  memory  materials.  Recall  from  [2]  that  this 
term  denotes  strain  histories  for  which  the  material  behaves  viscoelas- 
tically.  For  example,  in  a uniaxial  test  if  a permanent  memory  material 
is  taken  to  a strain  emax  and  then  brought  to  rest,  subsequently  the 
material  is  essentially  viscoelastic  as  long  as  strains  less  than  about 


For  notational  simplicity,  it  will  be  understood  that  in  formulas  such 
as  (2-1),  the  c's  may  be  time  dependent  kernals,  and  the  e's  may  be 
strain  histories.  Thus,  the  subsequent  analysis  will  be  valid  for 
viscoelastic  as  well  as  for  elastic  behavior. 
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75%  of  e are  applied.  Thus,  in  (2-2)  and  (2-3)  the  e..  are  not 
max  lj 

arbitrary,  but  must  be  in  a viscoelastic  range. (*) 

In  many  problems  of  practical  interest,  only  equations  (2-3)  are 
needed;  the  transient  formulas  of  Farris  [ 3 ] and  Quinlan  [ 4 ] provide 
more  information  than  necessary.  One  example  is  the  thermal  cycling 
problem:  a rocket  motor  is  subjected  to  a temperature  loading  and  then 

brought  to  rest.  The  thermal  loading  produces  damage;  therefore,  the 
material's  subsequent  mechanical  behavior  will  be  quite  different  from 
its  original  behavior.  It  is  important  to  predict  the  subsequent 
behavior  from  the  initial  modulus  and  the  thermal  loads;  that  is,  to 
determine  relations  (2-3).  A transient  description  of  the  material 
properties  during  the  thermal  loading  is  unnecessary. 

The  relations  (2-3)  can  be  obtained  as  special  cases  of  transient 
formulas  such  as  Farris'  or  Quinlan's.  However,  (2-3)  actually  states 
the  problem  more  generally  because  no  restriction  is  placed  on  how  the 
loads  L enter  the  equations.  Of  course,  this  generality  leads  to  an 
impractically  complex  theory,  and  simplifications  must  be  sought.  For 
example,  consider  the  simplest  special  case  of  (2-3):  the  material  is 

characterized  in  (2-1)  and  (2-2)  by  single  fourth  rank  modulus  tensors 
which  are  linearly  related  to  each  other.  Then  (2-3)  becomes 


If  (2-1)  is  used  to  describe  the  virgin  state,  the  viscoelastic  range 
may  reduce  to  infinitesimal  strains.  In  this  case,  c(0)  is  the  tangent 
modulus  at  zero  strain.  (2-1)  cannot  describe  the  hysteresis  effects 
in  the  virgin  state,  but  one  can  still  pose  the  problem  of  relating 
the  damaged  moduli  of  (2-2)  to  the  initial  tangent  modulus  through 
formula  (2-3). 
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(2-4) 


where  an  eighth  rank  tensor  function  D of  the  loads  L must  be  determined 
experimentally.  Even  allowing  for  the  index  symmetries  in  and  in 

(2-4)  requires  the  evaluation  of  an  excessive  number  of  constants. 
Furthermore,  the  dependence  of  D on  L is  still  unspecified. 

In  [ 1 ] Fitzgerald  suggested  that  damage  could  be  measured  by  a 
second  rank  tensor.  This  idea  is  motivated  by  observing  that  damage  could 
reduce  a material's  modulus  by  different  amounts  in  different  directions. 
Thus,  a damaged  initially  isotropic  material  could  well  be  anisotropic. 

A second  rank  tensor  is  the  simplest  quantity  which  can  model  this 
possibility.  Of  course,  more  complex  quantities  could  also  be  used: 
formula  (2-4)  illustrates  the  use  of  an  eighth  rank  tensor  to  model 
damage.  In  any  case,  in  general  damage  cannot  be  described  simply  by  a 
scalar . 


Analytically,  Fitzgerald's  suggestion  is  incorporated  into  equation 
(2-3)  by  replacing  L by  a second  rank  tensor  which  depends  on  L: 


,(D  . C(U  (c(0)  ) (*) 
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Equation  (2-5)  should  be  normalized  so  that  the  moduli  are  unchanged 
if  there  is  no  damage.  For  later  work,  it  will  be  convenient  to  let 
D^.  = 0 corresond  to  no  damage;  then 
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75%  of  e are  applied.  Thus,  in  (2-2)  and  (2-3)  the  e,,  are  not 
max  ij 

arbitrary,  but  must  be  in  a viscoelastic  range. (*) 

In  many  problems  of  practical  interest,  only  equations  (2-3)  are 
needed;  the  transient  formulas  of  Farris  [ 3 ] and  Quinlan  [ 4 ] provide 
more  information  than  necessary.  One  example  is  the  thermal  cycling 
problem:  a rocket  motor  is  subjected  to  a temperature  loading  and  then 

brought  to  rest.  The  thermal  loading  produces  damage;  therefore,  the 
material's  subsequent  mechanical  behavior  will  be  quite  different  from 
its  original  behavior.  It  is  important  to  predict  the  subsequent 
behavior  from  the  initial  modulus  and  the  thermal  loads;  that  is,  to 
determine  relations  (2-3).  A transient  description  of  the  material 
properties  during  the  thermal  loading  is  unnecessary. 

The  relations  (2-3)  can  be  obtained  as  special  cases  of  transient 
formulas  such  as  Farris'  or  Quinlan's.  However,  (2-3)  actually  states 
the  problem  more  generally  because  no  restriction  is  placed  on  how  the 
loads  L enter  the  equations.  Of  course,  this  generality  leads  to  an 
irapractically  complex  theory,  and  simplifications  must  be  sought.  For 
example,  consider  the  simplest  special  case  of  (2-3):  the  material  is 

characterized  in  (2-1)  and  (2-2)  by  single  fourth  rank  modulus  tensors 
which  are  linearly  related  to  each  other.  Then  (2-3)  becomes 


If  (2-1)  is  used  to  describe  the  virgin  state,  the  viscoelastic  range 
may  reduce  to  infinitesimal  strains.  In  this  case,  c^)  is  the  tangent 
modulus  at  zero  strain.  (2-1)  cannot  describe  the  hysteresis  effects 
in  the  virgin  state,  but  one  can  still  pose  the  problem  of  relating 
the  damaged  moduli  of  (2-2)  to  the  initial  tangent  modulus  through 
formula  (2-3). 
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(2-4) 
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where  an  eighth  rank  tensor  function  D of  the  loads  L must  be  determined 
experimentally.  Even  allowing  for  the  index  symmetries  in  and  in 

c(1),  (2-4)  requires  the  evaluation  of  an  excessive  number  of  constants. 
Furthermore,  the  dependence  of  D on  L is  still  unspecified. 

In  t 1 ] Fitzgerald  suggested  that  damage  could  be  measured  by  a 
second  rank  tensor.  This  idea  is  motivated  by  observing  that  damage  could 
reduce  a material's  modulus  by  different  amounts  in  different  directions. 
Thus,  a damaged  initially  isotropic  material  could  well  be  anisotropic. 

A second  rank  tensor  is  the  simplest  quantity  which  can  model  this 
possibility.  Of  course,  more  complex  quantities  could  also  be  used: 
formula  (2-4)  illustrates  the  use  of  an  eighth  rank  tensor  to  model 
damage.  In  any  case,  in  general  damage  cannot  be  described  simply  by  a 
scalar. 


Analytically,  Fitzgerald's  suggestion  is  incorporated  into  equation 
(2-3)  by  replacing  L by  a second  rank  tensor  D which  depends  on  L: 
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(2-5) 


Equation  (2-5)  should  be  normalized  so  that  the  moduli  are  unchanged 
if  there  is  no  damage.  For  later  work,  it  will  be  convenient  to  let 
D = 0 corresond  to  no  damage;  then 


(C(0)  0)  - c(0) 
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The  relation  between  D and  L is  still  unknown,  however,  the  general 
expansion  (2-4)  is  easily  worked  out.  The  details  are  carried  out  in 
sections  4 and  5,  and  some  consequences  of  (2-4)  are  discussed  there. 

It  should  be  emphasized  that  assuming  the  loads  L enter  (2-3)  through  a 
second  rank  tensor  is  a simplification  which  may  not  be  valid  under  all 
conditions.  The  character  of  the  loads  L may  be  important  in  determining 
the  validity  of  the  approximation.  This  issue  is  discussed  further  in 
section  3. 


2-3  Existing  experimental  evidence  does  not  suggest  how  L enters 
equation  (2-3)  when  L is  a completely  arbitrary  loading  history.  The 
reason  is  that  only  a special  class  of  loading  histories  is  usually 
produced  in  the  laboratory,  namely  those  histories  with  principal  axes  of 
stress  and  strain  constant  in  time.  Uniaxial  tests,  the  most  commonly 
performed  on  propellant,  are  obviously  of  this  type.  Assuming  that  damage 
is  produced  by  this  kind  of  loads  is  not  necessarily  too  restrictive 
since  practical  problems,  such  as  the  thermal  cycling  problem,  often 
involve  precisely  such  conditions. 


Keeping  in  mind  the  restriction  to  constant  principal  directions, 
experiments  do  suggest  that 

(a)  Damage  is  caused  by  tensile  strains. 

(b)  The  maximum  tensile  strain  is  the  most  important  factor  is 
causing  damage. 

Thus,  as  a first  approximation,  may  be  taken  as  the  maximum  tensile 
strain  during  the  loads  L.  A suitable  measure  of  this  strain  will  now 
be  developed.  Denote  the  strain  history  by  E(t),  and  let  E(t)  be  reduced 
to  principal  axes  by  the  time-independent  rotation  matrix  R.  Then 


E (t)  0 0 
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where  x,  y,  and  z denote  principal  axes.(*)  Let  E be  the  maximum  value 
of  E(t) . Then 
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is  a measure  of  maximum  tensile  strain  which  vanishes  if  the  principal 
strains  are  never  tensile . Alternatively,  damage  might  be  measured  by  some 
function  f(D).  If,  for  example,  tensile  strains  below  a certain  threshold 
Et  did  not  cause  damage,  one  might  choose  instead  of  (3-1), 


A , * 

E - E + E - E„ 
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Once  D is  related  to  L,  the  relations  (2-3)  between  the  initial  and  damaged 
moduli  are  obtained  by  using  the  results  of  sections  4 and  5. 

One  problem  which  does  not  involve  loading  histories  with  constant 
principal  directions  can  actually  be  handled  by  (2-3).  Suppose  that  the 


Recall  that  E = 6 . in  the  unstrained  state  and  that  principle  strains 

J ■*"  J 

E>1  are  tensile , while  principal  strains  E<1  are  compressive. 


t 1 


virgin  sample  with  moduli  is  subjected  to  damage  and  brought  to 

rest.  New  moduli  result.  Now  let  the  material  be  subjected  to 

another  damage  field  with  different  principal  axes  from  D^.  The 
resulting  moduli  c^  are  expressed  in  terms  of  the  c^  by  (2-3): 


, (D  n!\ 

(c,,  , D ) 

-*-J  * • • • • * 


To  express  the  c^2^  in  terms  of  the  c^\  use  (2-3)  to  express  the  c^^ 
in  terms  of  c^\  The  result  is 


(c 


(0) 

ijkl, 


Note  that  and  enter  this  equation  in  a complex  way.  There  need  not 
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be  any  one  second  rank  tensor  D such  that 
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This  example  indicates  the  complexities  possible  if  damage  is  produced  by 
general  loading  histories.  This  issue  is  not  pursued  any  further  here 
since  its  practical  significance  is  unclear. 


2-4  The  damage-dependent  stress-strain  relations  can  be  derived 
directly  from  equation  (2-5).  However,  it  is  simpler  to  proceed  by 
substituting  (2-5)  into  (2-2).  The  resulting  expression  for  stress  depends 
on  both  strain  and  damage: 


T = T(D,  E)  (4-1) 

Equation  (4-1)  is  really  equivalent  to  (2-5);  indeed,  if  D is  treated  as  a 
parameter  in  (4-1),  the  moduli  of  various  orders  can  be  derived,  and  the 
relations  (2-5)  recovered.  The  difference  between  the  two  equations  is 
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only  the  point  of  view.  In  (4-1),  only  D and  E are  variable,  whereas 
(2-5)  allows  the  moduli  c.^  to  vary  as  well.  Thus,  the  form  (2-5)  is 
required  only  when  the  moduli  are  changing;  typical  cases  are  superposing 
damage  fields  as  in  section  3,  or  developing  a transient  theory. 

The  easiest  approach  to  simplifying  (4-1)  is  based  on  polynomial 
approximation.  If  the  D„  and  the  do  not  exceed  a fixed  maximum, 
Weierstress'  theorem  implies  that  the  T may  be  taken  as  polynomials  in 
the  and  the  E and  the  difference  between  the  polynomial  form  and 
the  hypothetical  exact  form  may  be  made  arbitrarily  small.  Thus,  one 
can  set 


T.  , = E (E.  .)  P (D  )q 

ij  P,q  ljklmn  kl  mn 


(4-2) 


The  precise  range  on  p and  q in  (4-2)  is  not  specified,  but  it  is  assumed 

to  be  finite.  Equation  (4-2)  is  still  too  general  since  the  components 

T. . that  it  defines  may  not  transform  as  a second  rank  tensor.  To  assure 
ij 

the  proper  transformation  laws  T must  be  built  up  from  D and  E by  tensor 
operations , namely  by  tensor  products  and  by  contractions  with  constant 
tensors.  Invoking  these  principles,  (4-2)  becomes 


T = T(°0)  + T(*J}  D + T(?JJ  E + T(J?J  D E + ... 
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(4-3) 


where  is  a constant  2(i+j+l)t*1  rank  tensor.  Alternatively, 

(4-3)  may  be  rewritten  by  collecting  like  products  of  E.  Then 


T,,  = T,(0)  (D)  + T,]) , (D)  E.  . + T.^.\  (D)  E.,  E 
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where 
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It  should  be  emphasized  that  all  of  the  sums  in  (4-3),  (4-4),  and  (4-5) 
are  finite.  Of  course,  in  practice,  only  a very  small  number  of  terms 
would  be  taken. 

Equation  (4-4)  can  be  used  to  obtain  approximate  theories  of  various 
orders.  Consider  first  the  case  of  infinitesimal  strain.  Then  all 
products  involving  E can  be  ignored,  and  then  (4-4)  becomes,  writing 
e for  E: 


o = o(0)  (D)  + c(1)  (D)  e (4_6) 

In  (4-6),  o(0)  is  the  residual  stress  present  in  the  zero  strain 
configuration,  and  c^  (D)  is  the  damage  dependent  modulus  tensor. 

Next  consider  the  case  that  D is  small  and  the  strains  are  roughly 

an  order  of  magnitude  smaller  than  D.  This  case  would  apply  to  a test 

sequence  in  which  a material  is  subjected  to  strains  of  around  10% 

and  subsequently  tested  at  strains  of  one  or  two  percent.  Then  it  is 

2 3 

reasonable  to  retain  terms  of  order  D,  D , and  D and  DE  in  (4-6)  and  to 
discard  the  others;  thus,  (D)  is  at  most  cubic  in  D and  c^  is 

linear  in  D.  In  view  of  the  difficulty  in  measuring  the  residual 
stresses  accurately,  it  is  not  unreasonable  to  take  linear  in  D as 

well.  This  procedure  leads  to  the  simplest  theory  exhibiting  permanent 


stresses  and  modulus  reduction: 
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and  c 
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are  linear  in  D. 


(4-7) 


Many  test  conditions  closely  approximate  the  requirements  for  (4-7).  This 
fact,  together  with  its  simplicity,  makes  (4-7)  a practical  approximate 
theory. 


2.5  The  tensor  transformation  of  section  3 does  not  impose  any 
restriction  on  the  tensors  T^J)  appearing  in  (4-3).  The  form 
of  these  tensors  can  only  be  determined  from  considerations  of 
material  symmetry.  In  this  section,  the  stress-strain  relations  for  an 
initially  isotropic  material  will  be  developed.  For  such  a material, 
the  tensors  •••  of  (2-1)  are  isotropic;  therefore,  they  are  sums  of 

tensor  products  of  Kronecker  deltas.  It  seems  reasonable  that  the  T^^ 
will  also  be  isotropic  tensors.  If  so,  the  material  will  be  said  to 
exhibit  Isotropic  damaging. (*)  Thus,  if  strain  and  damage  are  rotated 
arbitrarily  while  the  material  is  fixed,  the  stress  is  rotated  the  same 
way.  On  the  other  hand,  it  is  not  impossible  for  an  initially  isotropic 
material  to  damage  anisotropically . For  example,  a given  tension  might 
produce  greater  degrees  of  weakening  in  some  directions  than  in  others. 
Whether  or  not  the  T^^  have  the  same  symmetry  properties  as  the  c^ 
can  only  be  determined  by  experiments. 

The  general  form  of  equation  (2-2)  for  an  initially  isotropic 
material  with  isotropic  damaging  is  found  by  replacing  the  T^^in 
ft-3)  by  arbitrary  products  of  Kronecker  deltas.  The  expression  that 
results  is  clearly  a sum  of  terms  of  the  form  tr(DmEn)  DP  Eq  . Rivlin 


has  shown  (Green  and  Adkins  [ 5 J , appendix)  that  this  sum  can  be  reduced 
to  the  form 


T 
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(5-2) 


where  the  are  polynomials  in  the  invariants 

2 3 2 3 

trE,  trE  , trE  , trD,  trD  , trD  , trED, 

(5-3) 

trE^,  trED2,  trE2D2 

Let  D be  fixed  in  (5-2)  . The  resulting  stress-strain  relation  has 
orthorhombic  symmetry  with  crystal  axes  coincident  with  the  principal 
axes  of  D.  To  verify  this  fact,  recall  that  (5-2)  satisfies  (5-1). 

Thus, 

RTTR  = T(D,  rter) 

T 

holds  whenever  R DR  = D.  Therefore,  (5-2)  with  D fixed  is  invariant 
under  symmetry  operations  which  leave  D invariant.  A reflection  through 
any  plane  formed  by  D's  principal  axes  leaves  D unchanged,  and  hence 
leaves  (5-2)  unchanged . (*)  Any  other  symmetry  operation  will  in  general 
change  the  form  of  D and  hence  change  the  form  of  (5-2).  Higher  degrees 
of  symmetry  are  possible  only  if  D is  more  symmetric.  Thus,  if  D has 
two  equal  principal  values,  it  is  invariant  under  all  rotations  in  the 
plane  formed  by  the  corresponding  principal  axes.  D,  and  therefore  (5-2) 


These  are  the  symmetries  characteristic  of  an  orthorhombic  material. 


are  then  transversely  isotropic.  Finally,  if  D has  three  equal  principal 
values,  it  is  invariant  under  all  rotations;  (5-2)  is  then  fully  isotropic. 

Thus,  imposing  a single  damage  field  transforms  an  isotropic 
material  into  an  orthorhombic  material.  Note,  however  that  the 
orthorhombic  symmetry  represented  by  (5-2)  is  of  a special  type.  For 
instance,  cubic  symmetry  cannot  occur  in  the  form  (5-2).  Interchanging 
the  principal  axes  of  D leaves  D invariant  only  if  D has  three  equal 
principal  values;  in  this  case,  (5-2)  is  fully  isotropic.  This 
conclusion  also  results  from  a count  of  arbitrary  constants.  The  general 
linear  orthorhombic  material  is  characterized  by  nine  constants.  In 
(5-2),  however,  the  arbitrary  constants  in  a linear  theory  are  the  two 
initial  isotropic  moduli,  and  the  three  principal  values  of  D,  for  a 
total  of  only  five  constants. 

2.6  Equation  (5-2)  contains  eight  arbitrary  polynomials  f^  ...  fg. 
Therefore,  theories  containing  any  number  of  free  constants  are 
available  for  material  characterization  in  the  damaged  regime.  Once  a 
definite  theory  has  been  selected,  its  unknown  constants  are 
evaluated  by  curve  fitting  to  experimental  data  on  laboratory  tests  such  as 
uniaxial  tension.  This  procedure  is  no  different  in  principle  from 
that  used  in  finite  elasticity  (see  for  example,  [5],  ch.  11  ).  Just 
as  in  finite  elasticity,  care  must  be  taken  to  select  a consistent 
theory.  For  example,  if  a quadratic  term  in  E appears  in  the  theory, 
then  all  possible  quadratic  terms  in  E must  appear  as  well.  Tests  may 
reveal  that  some  of  these  additional  terms  have  zero  coefficients,  but 
this  situation  must  be  established  by  experiments  and  should  not  be 
assumed  a priori. 
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The  use  of  (5-2)  in  material  characterization  is  simplified  somewhat 
since  its  terms  correspond  to  a few  easily  identified  qualitative  effects 
in  the  damaged  state.  Accordingly,  (5-2)  will  now  be  analyzed  in 
more  detail  from  this  point  of  view.  In  what  follows,  p(D)  and  q(E) 
will  denote  polynomials  in  the  scalar  invariants  of  D and  E,  not 
necessarily  the  same  ones  each  time  they  appear.  Also,  it  will  be 
convenient  to  write  the  stress-strain  relations  in  terms  of  a strain 
measure,  for  example  E - I,  which  vanishes  in  the  undamaged  unstrained 
state. 

2 3 ? 3 

(a)  Terms  independent  of  E:  p(trD,  trD  , trD  )I,  p(trD,  trh-,  trD  )D, 

p(trD,  trD2,  trD^)D2 

Let  the  damaged  material  be  brought  to  its  original  unstrained  state. 

These  terms  then  appear  in  (5-2)  as  stresses  necessary  to  do  this. 
Measurements  of  these  stresses  would  thus  be  used  to  evaluate  coefficients 
which  correspond  to  the  terms  ( a ) . Since  the  stress  T is  measured 
from  the  initial  rest  configuration,  p I is  a purely  hydrostatic  stress, 
while  the  other  two  terms  contribute  stresses  which  depend  on  the 
principal  directions  of  D.  Evaluating  the  polynomials  p(D)  will,  of 
course,  require  experiments  with  different  degrees  of  damage  D. 

Terms  ( a ) may  be  regarded  alternatively  as  measures  of  permanent 
set : this  is  the  strain  present  in  the  damaged  configuration  under 

zero  stress. (*)  To  evaluate  this  strain,  let  T = 0 in  (5-2)  and  solve  for 


More  precisely,  "strain"  means  strain  measured  with  respect  to  the 
initial  rest  state,  and  by  "zero  stress,"  zero  surface  traction  is 
understood. 
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E.  Because  of  terms  ( a ),  E * 0 is  not  in  general  the  solution. 
Theoretically,  this  interpretation  makes  it  possible  to  evaluate  coeffi- 
cients corresponding  to  ( a ) by  measuring  permanent  set,  rather  than 
the  stresses  needed  to  remove  it.  However,  this  process  requires  the 
solution  of  (5-2)  with  T = 0;  in  effect,  (5-2)  must  be  solved  for  E 
in  terms  of  T and  D.  This  solution  will  be  practical  only  in  very 
simple  cases. 

(b)  p(D)q(E)  E,  p (D)q (E)  E2 

D enters  these  terms  through  its  scalar  invariants  only.  The  principal 
directions  of  D therefore  play  no  role.  If  just  terms  ( b)  appear  in 
(5-2),  an  initially  isotropic  material  remains  isotropic  after  damage. 
Thus,  if  damage  was  caused,  for  example,  by  a uniaxial  tension,  subsequent 
uniaxial  tests  would  give  identical  results  in  all  directions  of  test. 

The  direction  of  the  damage  field  would  not  be  distinguished.  This 
situation  can  be  called  isotropic  damage . 

Qualitatively,  a number  of  different  effects  can  be  modeled  with 

these  terms.  Let  the  initial  stress-strain  relation  be 

2 

T = T(E)  = qQ(E)I  + q^(E)E  + q2(E)E  . If  subsequent  to  damage, 

T = [1  + p(D)]T(E),  all  stress-strain  curves  are  reduced  proportionally; 
that  is,  the  damaged  curve  for  any  test  is  a constant  multiple  of  the 
undamaged  curve.  Furthermore,  the  constant  of  proportionality  is  the 
same  for  all  tests,  and  depends  only  on  the  degree  of  damage.  This  is 
the  simplest  theory  of  isotropic  damage. 

A more  general  theory  results  from  setting 
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T = qo(E)[l  + Po(D)]I  + q2(E)[l  + p2(D)]E  + q2(E)[l  + p2 (D) ] E2 

In  this  case,  the  damaged  and  undamaged  curves  for  the  same  test  need 
not  be  proportional.  Finally,  if  arbitrary  terms  ( b ) appear, 
different  cross-effects  can  be  present  in  the  damaged  and  undamaged 
states . 

(c)  q (det  E)p (D) I , q(det  E)p(D)D,  q(det  E)p(D)D2 

First,  note  that  since  det  E is  a polynomial  invariant  of  E,  it  can  be 

2 3 

written  as  a polynomial  in  tr  E,  tr  E , and  tr  E . When  q(l)  = 0,  these 
terms  have  the  following  interpretation.  Let  the  damaged  material  be 
returned  to  its  initial  unstrained  state.  Any  stresses  needed  to  do 
this  are  ignored  in  what  follows.  These  stresses  would  be  covered  under 
( a ) in  any  case.  Now  subject  the  damaged  sample  and  an  undamaged 
sample  to  identical  volume-preserving  deformations . (*)  Since  then 
det  E = 1,  and  q(l)  = 0,  when  only  terras  ( c ) appear  in  (5,2),  the 
stresses  in  the  undamaged  and  damaged  bodies  are  the  same.  Thus,  there 
is  no  difference  in  mechanical  response  to  a volume  preserving 
deformation. 

When  a deformation  which  is  not  volume  preserving  is  applied,  terms 
( c ) represent  additional  stresses  required  to  impose  the  deformation 
in  the  damaged  state.  The  first  term  in  ( c ) is  distinguished  from  the 
other  two  by  its  behavior  under  volumetric  strain  tests.  If  such  a test 
requires  just  hydrostatic  stress,  then  only  the  first  term  is  present  in 


* 


"Volume  preserving"  means  volume-preserving  with  respect  to  the 
undamaged,  unstrained  state. 


the  stress-strain  relation.  If  the  stress  is  not  hydrostatic,  then  the 
other  two  terms  are  present.  Note  that  in  this  case,  the  principal  axes 
of  stress  must  coincide  with  the  principal  axes  of  D. 

These  terms  are  not  appropriate  for  materials  which  are 
incompressible  in  both  the  damaged  and  undamaged  state,  since  then 
det  E is  constant. (*) 

When  q(l)  ^ 0 or  when  det  E is  replaced  by  other  invariants  of  E, 
then  terms  ( c ) have  analogous,  but  physically  somewhat  implausible 
interpretations.  Thus,  if  q(a)  = 0,  the  damaged  and  undamaged  materials 
respond  identically  to  deformations  such  that  det  E = a.  If  also, 
det  E is  replaced  by  some  invariants  I^(E),  and  q(a^,  ...  a^)  = 0,  then 
the  mechanical  response  is  identical  when  I^CE)  = a^.  If  q is  never 
zero,  terms  ( c ) simply  represent  additional  stresses  needed  to  impose 
certain  deformations  on  the  damaged  state;  these  stresses  depend  only  on 
certain  invariants  of  the  deformation.  These  terms  would  seem  to  be  of 
most  interest  in  the  case  discussed  at  first. 

Finally,  to  obtain  the  corresponding  terms  in  a linear  theory  with 
infinitesimal  strains,  det  E is  replaced  by  tr  E,  and  necessarily 
q(x)  ■ Ax,  otherwise  the  theory  is  not  linear.  Then  terms  ( c ) have 
the  additional  interpretation  that  the  damaged  and  undamaged  shear  modulus 
are  the  same. 


* 


Of  course,  the  constant  need  not  be  the  same  in  both  the  undamaged  and 
damaged  states  since  damage  may  cause  an  increase  in  volume. 


(d)  Terms  containing  products  of  E and  D and  their  invariants. 


These  terms  cause  anisotropic  behavior  in  the  damaged  state.  For 

example,  let  = D,  all  other  D^  * 0;  this  damage  field  would  be 

caused  by  extension  in  the  x direction.  Now  let  E represent  an  arbitrary 

uniaxial  tension.  Then  trDE  has  a maximum  for  tension  in  the  x direction 

and  will  have  minimums  for  tension  in  the  y and  z directions.  This  fact 

makes  it  easy  to  determine  whether  or  not  tr  DE  should  appear  in  the 

damaged  stress-strain  relation.  The  other  joint  invariants  of  E and  D: 

2 2 2 2 

tr  E D,  tr  ED  , and  tr  E D , have  similar  properties  and  can  be  used 
for  curve  fitting. 


A simple  example  of  ( ) is  T = (1  + a tr  DE)  T(E)  where,  as  in 

( c ) , T = T(E)  is  the  undamaged  stress-strain  relation.  In  this 
theory,  the  stress-strain  curves  in  the  damaged  and  undamaged  material 
are  proportional  for  any  test  procedure.  However,  the  constant  of 
proportionality  depends  on  the  relative  orientation  of  E and  D.  As  in 
( c ),  more  complex  formulas  of  this  type  can  be  constructed  quite 
easily.  Thus,  in 


T = qQ(E)[l  + Pq(d)  ] + q1(E)[l  + p^DjJE  + q2(E)[l  + P2(D)]E2 


the  damaged  and  undamaged  curves  are  not  proportional,  and  the  degree 

* 

to  which  they  differ  depends  on  the  orientation  of  D and  E. 

The  terms  ED  + DE,  etc.  act  similarly  when  only  tension  tests  are 
involved.  They  can  be  distinguished  from  the  terms  involving  joint 
invariants  only  by  their  behavior  under  volumetric  straining  of  the 
damaged  material.  Let  E = a I.  Then  the  joint  invariants  are  multiples 
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— 


of  the  invariants  of  D.  When  only  these  terms  occur  in  (5-2),  the 
principal  directions  of  D have  no  effect  on  the  damaged  material,  and  a 
volumetric  strain  requires  only  hydrostatic  stress.  But  if  the  terms 
ED  + DE  etc.  are  present,  the  volumetric  strain  will  require  a non- 
hydrostatic stress. 


For  incompressible  materials,  the  volumetric  strain  is  inadmissable. 
Then  the  two  terms  can  be  distinguished  by  biaxial  damage  fields.  Let 


D 

E 
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yy 

yy 
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When  only  ED  + DE  etc.  occur  in  (5-2),  the  stresses  in  the  x and  y 
directions  are  independent,  whereas  the  joint  invariants  introduce 
coupling . 


( a ) through  ( d ) exhaust  all  of  the  terms  which  can  occur  in 
(5-2).  Summarizing  their  properties,  terms  ( a ) are  required  whenever 
permanent  set  is  observed.  When  damage  is  isotropic,  terms  ( b ) are 
sufficient  for  a characterization.  Isotropic  damage  is  easily  detected 
from  tension  tests  in  different  directions.  The  description  of 
anisotropic  damage  requires  terms  ( d ).  Finally,  ( c ) describes 
specialized  situations  when  the  damaged  and  undamaged  materials  respond 
identically  to  certain  types  of  deformation. 


2.7  The  derivation  of  stress-strain  relations  in  sections  4 and  5 
follows  Cauchy's  approach  of  relating  the  stress  and  strain  tensors 
directly.  The  alternative  approach  of  Green,  in  which  the  relations  are 
obtained  through  a strain  energy  function,  will  be  outlined  in  this 
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section.  This  formulation  is  especially  useful  in  a coupled  thermo-mech- 
anical analysis  since  then  the  strain  energy  function  need  only  be 
replaced  by  the  Helmholtz  function. 


Let  the  initial  strain  energy  function  be  W0(E^).(*)  Then  the 
initial  stress-strain  relation  is 


3W 

T = n — 2 
ij  po3E 


ij 


After  a damage-producing  load  sequence  L,  the  strain  energy  is  W(E^,  L)  . 
As  in  section  4,  Fitzgerald's  suggestion  that  damage  be  measured  by  a 
second  rank  tcnso’"  D leads  to  the  strain  energy  W(E^j  , D^).  The 
relation  between  D and  L has  already  been  discussed  in  section  3.  As 
in  section  4,  when  the  components  E and  D^  do  not  exceed  a fixed 
maximum,  Weierstrass ' theorem  may  be  invoked  to  approximate  W by  a 
polynomial  in  the  E^  and  . No  further  restriction  on  W is  possible 
without  bringing  in  material  symmetry  properties.  If  a material  exhibits 
isotropic  damaging  (see  section  5 ) , W is  then  an  isotropic  function  of 
E and  D.  Then  Rivlin's  reduction  (see  [5  J,  p.  317)  leads  to  the 
conclusion  that  W is  a polynomial  in  the  basic  invariants 


Ix  = tr  E 
J1  = tr  D 


I9  = tr  E 

z 2 
« tr  D" 


tr  DE  K„ 


2 

tr  E D 


I,  = tr  EJ 

J 3 

J.  = tr  D 

J 2 
K„  = tr  ED 


K,  = 


2 2 

tr  D E 


(7-1) 


WQ  may  be  the  potential  for  the  initial  tangent  modulus  at  zero  strain. 
See  section  3. 
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The  damage  dependent  stress-strain  relation  is 


_ 3W  ^1  aw  aK4 

T1J  ' »VElj  + "•  5K4  aEtJ 

The  partial  derivatives  of  the  invariants  with  respect  to  are  as 
follows : 


Substituting  into  (7.2), 


ij 


3W 

31 


1 


ij 


+ 2 


3W 

31, 


i j 


+ 3 


3W 

31. 


(E2) 


ij 


oW 

3k, 


ij 


3W_ 

3k. 


(ED  + DE) 


ij 


3W 

3k. 


(D2) 


ij 


ls_ 

3k, 


(ED 


D2E) 


ij 


As  expected,  this  is  a special  case  of  formula  (5.2). 

2.8  The  stress-strain  relations  have  been  developed  using  the 
undeformed  body  as  the  reference  configuration;  this  is  the  Lagrangian 
formulation. (*)  In  this  section,  the  alternative  formulation  using  the 
deformed  configuration  as  the  reference  (the  Euler ian  formulation)  is 
developed . 


In  finite  elasticity,  the  choice  between  the  Eulerian  and  Lagrangian 
viewpoints  is  basically  arbitrary.  In  the  theory  of  permanent  memory 
materials,  however,  the  Lagrangian  viewpoint  has  the  definite  advantage 
that  the  damage  tensor  D,  defined  in  the  undeformed  undamaged  body,  is 
constant.  When  the  deformed  body  is  the  reference  configuration,  the 
damage  measure  depends  on  the  current  configuration.  For  this  reason, 
the  Lagrangian  formulation  will  be  taken  as  the  basic  formulation,  and 
the  Eulerian  equations  will  be  derived  formally  from  it. 


Recall  that  the  basic  equations  for  isotropic  damaging  were: 


T = f , I + f„D  + f.E  + f,D2  + fcE2  + f c (ED  + DE) 
1 2 3 4 5 6 


(8.1) 


22  2 2 2222 
+ fy(ED  + D E)  + fg(E  D + DE  ) + f^(E  D + D E ) 


More  precisely,  the  und amaged  undeformed  configuration  is  the  reference. 
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•X-  f 


where  the  f are  polynomials  in  the  invariants  (5-3).  T is  the  symmetric 

Kirchoff  stress,  and  E is  the  Green  strain  defined  in  Cartesian  coordinates 
2 

by  ds  = dX.dX..  The  Eulerian  equations  use  the  Cauchy  stress  t and 

2 

the  strain  e defined  by  dS  = e,. dx.dx.  . Before  transforming  to 

ij  i J 

Eulerian  variables,  it  is  important  to  decide  whether  (8-1)  contains 
contravariant,  covariant,  or  mixed  tensors.  Ordinarily,  this  distinction 
is  not  important  since  the  tensor  types  coincide  when  only  Cartesian 
coordinates  are  used.  However,  making  a transition  from  Lagrangian  to 
Eulerian  variables  generally  involves  transformation  to  an  oblique 
curvilinear  coordinate  system,  even  if  Cartesian  coordinates  are  used  in 
the  reference  configuration.  Thus,  the  distinction  is  important  in  this 
context . 

The  simplest  formulas  result  if  (8-1)  is  written  in  terms  of  mixed 

tensors.  Since  a self-contained  comprehensible  discussion  of  tensor 

analysis  is  impractical  we  list  the  important  properties  of  mixed 

tensors  and  refer  to  [ 6 ] for  details.  Mixed  components  of  a covariant 

tensor  A.,  are  defined  by 
ij 


*1  • ^ 


(8-2) 


where  dS2  = G . . dX1  dX-’ , and  Gik  G.  = 5*  . If  x = x (X.)  defines  a 
ij  kj  j i i : 

change  of  variables,  A*  transforms  as 


3X1  . 

j 9\  9x.  A1 


(8-3) 


Denote  by  X^  and  x^  Cartesian  coordinates  in  the  undeformed  and  deformed 
bodies  respectively.  In  terms  of  the  deformation  gradients 
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J M,  'j  S*k 


(8-  3)  becomes 


A*'  - f*  f1 
J k j 1 


or  in  matrix  notation 


A'  = FAf 


Since,  in  view  of  the  chain  rule  F*  = ff  F.  = 6^, 

k j k j j’ 

f = F * . Thus , the  mixed  components  transform  like  a matrix  under  a 
change  of  variables . (*)  Thus,  mixed  components  are  particularly  suitable 
for  matrix  equations  such  as  (8-1).  Indeed,  expressions  such  as 


(E2)*  - E*  Ej,  (ED)*  = E*  eJ,  trE  = E^ 


are  valid  tensor  operations,  while  corresponding  covariant  or  contr- 
variant  expressions  are  more  complex. 

The  mixed  Cauchy  stress  tensor  t*  is  related  to  T*  by 

Jt*  = F,1  f1  T^  (JT  = FTf ) 

1 k j 1 


where  J = det  F.  Define  the  Eulerian  damage  measure  d^  by 

d*  = F*  f J D*  (d  = FDf ) 

As  expected,  the  Eulerian  damage  depends  on  the  current  configuration. 
Next,  note  that 


This  fact  makes  it  particularly  convenient  to  define  the  mixed  damage 
tensor  D*-  in  terms  of  the  Ej . Recall  from  section  3 that  D was  defined 
as  a matrix  function  of  the  strain  history  E. 


r * 

r i 

'j.  : i 


T T -1 

FEf  = FF  Ff  = FF  = e 

Therefore,  if  p(D,E)  denotes  any  polynomial  in  D and  E appearing  in 

(8-1), 

Fp(D,  E)f  = p(d,  e_1)  (8-6) 

It  follows  from  this  equation  that  the  joint  invariants  transform 
according  to 


tr  p (D , E)  = tr  Fp(D,  E)f  = tr  p(d,  e ^) 


(8-7) 


Substituting  (8-5)  - (8-7)  into  (8-1), 


T. i ri  , ,i  , +i  , .,2,1  , , -2,  i 

Jtj  = 81  6j  + g2dj  + 836j  + 84(d  >j  + 85(e  >j 


(8-8) 


, -1J  ^ , -l.i  . , -1 , 2 2 -1  i . , -2,  ^ , -2. i , -2,2  L ,2  -2. 

+ g,(e  d + de  ).  + g7(e  d +d  e ) + g (e  d + de  ),  + gQ(e  d + d e ) 

o i / jo  j y 


where  the  g^  result  from  replacing  E by  e and  D by  d in  the  f^. 

Eulerian  stress-strain  relations  with  covariant  or  contravariant 
components  are  obtained  from  (8-8)  by  raising  or  lowering  indices  with 
the  metric  tensors  g^  and  g • But  it  is  important  to  bear  in  mind  that 
since  (8-8)  is  a transformation  of  (8-1),  these  are  the  metrics  induced 
in  the  deformed  body  by  the  deformation  (X ^ ) when  the  undeformed  body  has 

the  usual  Euclidean  metric.  Thus, 


S1J-  * F^-Ce-Vhs  - I 

1<  a<  3 J 1<  a<  3 J J 


This  situation  causes  the  contravariant  and  covariant  Eulerian  equations 
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'I 


to  be  more  complicated  than  (8-1).  Since  these  expressions  will  not  be 


used  in  what  follows,  rather  than  working  them  out  in  detail,  we  just  cite 

2 i 

one  typical  case:  suppose  that  one  wants  to  raise  indices  in  (d  )^. 


Multiplying  by  g 


U 


.i  ,k  lj  .1  ,kj  .il  ,kj  ..ij 

\ di  8 = dk  d = d 8lk  d = (ded) 

2 

Thus,  in  the  contravariant  form  of  (8-8),  d is  replaced  by  ded,  and 
2 

tr  d by  tr  ded. 

2.9  So  far,  the  finite  stress-strain  relations  have  been  expressed 
in  terms  of  the  original  undeformed  state,  and  in  terms  of  the  current 
deformed  state.  With  permanent  memory  materials,  it  is  also  possible 
to  state  the  relations  in  terms  of  the  damaged  traction-free  state. 

This  alternative  does  not  lead  to  a universal  form  of  the  stress- 
strain  relations,  as  do  the  Eulerian  and  Lagrangian  equations.  The 
reason  is  that  the  damaged  traction  free  state  depends  on  the  geometry 
of  the  body  and  on  the  damage  field.  In  view  of  this  fact,  we  will  only 
outline  the  procedure  for  obtaining  these  equations. 


First,  it  is  necessary  to  determine  the  damaged  traction  free 
state.  This  requires  solving  the  equilibrium  equations 


strain-deformation  relations 


E 


ij 


(9-1) 


(9-2) 
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k*  *■  k 


and  stress-strain  relations  (5-2)  subject  to  the  boundary  condition  of 

zero  surface  traction.  Because  of  the  terms  independent  of  E in  (5-2), 

the  solution  need  not  be  E = I,  T * 0.  Denote  the  solution  for  the 

deformations  by  x^  and  the  corresponding  deformation  gradient  by  F.  The 

problem  is  to  reformulate  (5-2)  in  terms  of  the  new  configuration  defined 

by  the  x^.  This  problem  is  analogous  to  the  transition  to  Eulerian 

variables  of  section  8;  the  only  difference  is  that  F is  a constant 

deformation  gradient  field,  so  that  the  strains  transform  differently. 

Denote  the  mixed  components  of  stress,  strain,  and  damage  in  the  traction 

free  state  by  T*,  E^  , and  D*  . Then 
•^1  ~3 


T1  = F*  T^ 
—1  k 1 1 


' Fk  fl  E1 
= f j Dk 

Dj  - Fk  fx  Dx 


Substituting  into  (5-2), 


4 ■ T1 4 


Recall  that  coordinates  have  been  so  chosen  that  T = 0 on  the  boundary 

- n 


of  the  reference  configuration  when  IJ  = I. 


To  raise  or  lower  indices  in  (9-4),  the  metric  tensors  G„  and  G 
must  be  used.  As  in  section  8,  these  are  the  metrics  induced  in  the 
damaged  state  by  the  deformation  when  the  undamaged  state  has  the 
Euclidean  metric  dS^  = dX^  + dY^  + dZ^.  Thus,  (9-4)  is  based  on  a 
coordinate  system  which  is  in  general  oblique  and  curvilinear. 


50 


It  is  also  possible  to  obtain  analogs  of  (9-4)  in  which  G*^  and  G. 
have  a simpler  form.  The  optimal  form  for  the  metric  will  depend  on  the 
geometry  of  the  damaged  state;  suppose,  for  example,  that  it  is  desired 
to  use  the  Euclidean  metric.  Then  the  metric  in  the  undamaged  state 
must  be  altered.  Let  G denote  this  metric.  G and  G are  related  by 


Therefore,  if 


kl 


1<  a<  3 


fk  fl 
a a 


then,  as  required. 


6 


ij 


Of  course,  formally  the  stress-strain  relation  so  obtained  are 
identical  to  (9-4);  only  the  coordinate  system  is  different.  Although 
indices  can  be  raised  and  lowered  using  5 and  <5^  so  that  the 
distinction  between  the  tensor  types  disappears,  equation  (5-2)  must  now 
be  rewritten  to  reflect  the  fact  that  the  metric  in  the  undamaged  state 


is  not  Euclidean.  For  example,  E written  in  terms  of  covariant 


components  is  E^  E^  G‘ 


kl 


2.10  For  some  applications,  the  finite  deformation  theory 
developed  so  far  may  be  unnecessarily  complex;  this  section  develops 
simple  linear  equations  for  permanent  memory  materials  in  a viscoelastic 


range . 


The  most  systematic  approach  to  deriving  linear  theories  starts  with 


the  exact  equations  and  introduces  suitable  approximations,  such  as 
infinitesimal  strains,  which  lead  to  linearity.  Although  this  approach 
will  not  be  followed  here,  it  is  worthwhile  to  understand  what  would  be 
involved  in  carrying  it  out.  The  first  linearization  of  equation  (5-2) 
would  not  be  to  introduce  inf initessimal  strains  directly.  Instead, 
one  would  start  with  the  theory  of  section  9 and  assume  infinitesimal 
strains  about  the  damaged  traction-free  configuration.  This  would  be 
an  example  of  infinitesimal  strain  superposed  on  finite  strain  (see  [ 5 ], 
§3.1  ) since  the  strains  needed  to  deform  the  original  stress-free  state 
to  the  damaged  state  might  be  finite.  Although  this  theory  would  be 
linear  in  the  strain  measures,  using  it  would  still  require  solution  of 
the  nonlinear  boundary  value  problem  which  determines  the  damaged 
traction-free  configuration. 


A second  degree  of  approximation  is  to  assume  that  the  strains  in 
the  damaged  traction  free  configuration  are  inf initessimal . This  is  the 
theory  which  will  be  developed  in  this  section.  This  theory  is 

(a)  linear  in  D and  in  the  infinitesimal  strains  E_^ 
and  more  importantly, 

(b)  the  total  strain  in  the  damaged  body  is  the  sum  of  the  strain 
needed  to  bring  the  original  rest  state  to  the  damaged  traction- 
free  configuration  and  the  strain  needed  to  bring  the  damaged 
traction  free  configuration  to  the  present  configuration. 


This  type  of  theory  is  appropriate  when  damage  is  caused  by  strains 
of  the  order  of  10%  to  20%,  and  subsequent  testing  involves  strains  of 
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a few  percent. (*)  Thus,  since  the  permanent  sets  are  rarely  large,  the  total 
strains  are  certainly  lesr  than  257.  and  infinitesimal  strain  theory  is 
adequate . 

Damage  is  defined  as  before  in  terns  of  the  maximum  tensile  strain: 
d = - max  [e(t)  + |c(t)|]/2 

where  |e(t)|  is  the  matrix  absolute  value  function.  The  stress-strain 
relations  are 


o = o(d,  e) 

The  ideas  of  isotropic  damaging  developed  in  section  2,  and  the  use  of 
Rivlin's  expansions  as  in  section  4 are  still  appropriate.  However,  in 
using  the  expansion,  all  terms  which  are  nonlinear  in  e are  dropped.  The 
question  then  arises  of  which  of  the  other  terms  ought  to  be  retained. 
Because  the  characterization  applies  to  a viscoelastic  range,  d must  be 
several  times  larger  than  e;  thus,  it  is  reasonable  to  retain  some 
products  of  d and  e and  even  perhaps  some  powers  of  d even  though  d is 
much  smaller  than  one.  Also,  it  is  desirable  to  model  anisotropy  in 
the  damaged  material.  On  the  other  hand,  the  experimental  data  do  not 
really  justify  a very  complex  theory.  In  view  of  these  considerations, 
the  following  equation  is  proposed  for  an  initially  isotropic  material 
with  Lam4  constants  X and  y: 


The  subsequent  strains  must  be  smaller  to  Insure  that  the  material 
remains  in  a viscoelastic  range. 
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The  remaining  constants  In  (10-2)  model  modulus  changes  a^  and  a^ 
model  Isotropic  damage  in  which  an  initially  isotropic  material  remains 
isotropic  after  damage,  but  with  possibly  reduced  moduli.  Note  that  a ^ 
models  a change  in  bulk  response  only. 


a22  corresponds  to  terms  (c)  , section  6 . When  tr  e = 0,  that  is, 
for  a volume  preserving  deformation,  this  term  does  not  contribute  to 
the  damaged  stress-strain  relation.  Thus,  it  models  an  anisotropic 
change  in  bulk  response  only;  the  shear  moduli  are  unchanged. 


The  other  constants,  a^  and  a^,  result  in  orthorhombic  behavior 
in  the  damaged  state.  These  terms  modify  both  bulk  and  shear  properties. 

Alternatively,  (10-1)  and  (10-2)  can  be  derived  from  a damage-dependent 

strain  energy  function  W = W(d,  e)  as  in  section  7 . To  derive  the 

type  of  theory  developed  here,  W should  be  at  most  quadratic  in  e and 

12  2 

linear  in  d.  Also,  W should  reduce  to  y A (tre)  + ptr(e  ) when  d = 0. 
Thus , 


W = y A (tre)2  (1  + a12  trd)  + ptre2  (1  + a31  trd) 


+ a21  tr(ed)  + a^  tre  d = a ^ (tre)  (trd)  + a22  tr(ed)  tre 


Substituting  into  a 


3w 


id  ^ 


leads  to 


o = [a^^  trd  + A (1  + a12  trd)  tre  + a22  tr(ed)]  I 


+ (a21  + a22  tre)  d + 2u(l  + a31  trd)  e + a^  (ed  + de) 


(10-3) 
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This  is  the  same  as  equation  (10-1)  and  (10-2),  except  that  = a ^ . 

To  rewrite  Q.Ci-3)  in  the  form  a..  = + C.  e,  , choose  axes  so 

ij  ij  ijkl  kl 

that  d is  in  diagonal  form.  Then  the  elastic  constants  of  the  damaged 
material  are: 


Clill  'i(1  + an  trd>  + 2H(1  + a,,  trd)  + (a,  , + ao„  + 2aA1  ) d 


12 


31 


13  22  41  11 


C2222  “ X^1  + a io  trd)  + 2y^1  + aii  trd)  + (a,o  + a„  + 2a,.,)  d 


12 


31 


13  22  41  22 


C3333  = X^  + aio  trd)  + + an  trd)  + (an  + aoo  + 2a,.,)  d 


12 


31 


13  22  41'  33 


C1122  = X(1  + a!2  trd)  + a22(dll  + d22> 


C1133  X(1  + al2  trd)  + a22(dll  + d33) 


C2233  = X(1  + al2  trd)  + a22(d22  + d33} 


C1212  " 2M(1  + a31  trd)  + 2a41(dll  + d22) 


C1313  2m(1  + a31  trd)  + 2a41  (dll  + d33) 


C2323  = 2u(1  + a31  trd)  + 2a41  (d22  + d33) 


These  nine  constants  are  characteristic  of  an  orthorhombic  material. 

The  crystal  axes  coincide  with  the  principal  axes  of  d.  However,  (10-3) 
does  not  represent  a general  orthorhombic  material.  Indeed,  as  in 
section  5 , (10-3)  cannot  have  cubic  symmetry  since  cmi  = c2222  = 

C3 3 J 3 ^orces  dn  = d22  = d33*  leads  to  isotropy. 
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List  of  Symbols  for  Chapter  III 


c 

Green's  tensor 

a 

Duration  of  a process 

p 

y-.o 

Viscosity  metric 

Maximum  metric 

V*0 

Weighted  Lebesgue  metric 

ar(-.*) 

Relaxation  metric 

bg(-  »■) 

Configuration  metric 

y(-,-) 

Metric  on  the  state  space 

a (•  ,*) 

History  metric 

Relaxation  History  metric 

aD('  .■) 

Damage  history  metric 

^ q( * ’ ’ ) 

Configuration  history  met 

D Damage  function 

D Damage  tensor 

D*  Relative  damage  history 


Damage  resulting  from  a constant  strain  E 
Damage  space 
Damage  functional 
Function;  influence  function 
G Configuration 

G Configuration  mapping 

G Tensor  valued  kernel 


G* 


I 


I 

Y 


A 

A 

o 

n 

P° 

n 

nG 

R 

P 

P 

s 

S 

S(H) 

8 

o 

Z 


ZG 


Stress-free  configuration 

Freeze  for  time  t at  configuration  G 

Configuration  space 

Strain  invariant 

Volume  dilation 

Octahedral  shear  strain 

Absolute  value 

pth  order  Lebesgue  average 

State  relaxation  mapping 

Virgin  state 

Unit  normal 

Set  of  processes  starting  at  state  a and  ending  in 
configuration  G 

Set  of  deformation  processes 

Set  of  processes  that  start  at  configuration  G 

Rotation  tensor;  strain  rate 

Evolution  function 

History  evolution  function 

Stress  value  of  the  sample;  dummy  variable 

Stress  mapping 

Stress  corresponding  to  a history  H 
Response  functional 
Relaxed  state 
State  space 
G section  of  z 
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CHAPTER  III.  GENERAL  THEORY 


3 . 1 Introduction 

The  hysteresis  erfect  associated  with  Permanent  Memory  materials 
can  be  visualized  as  a weakening  or  damaging  of  the  material  as 
demonstrated  in  Sul i joadikusumo  and  Fitzgerald  [ 1 ].  Quinlan  and 
Fitzgerald  [2.31  have  presented  this  concept  in  a constitutive 
equation  and  incorporated  it  into  a characterization  applicable  to 
highly  filled  polymeric  materials.  They  did  not,  however,  define 
an  order  relation  on  damage.  This  was  accomplished  by  Fitzgerald  [4  ]. 

This  present  paper  expands  upon  and  refines  the  work  of  Fitzgerald 
[ 4].  The  structure  developed  by  Noll  [ 5 ] and  Fitzgerald  [ 4 ] will 
be  used  herein. 

The  topology  on  the  state  space  £ is  first  developed,  followed 

by  a discussion  on  some  metrics  possible  for  this  space.  It  is 

shown  how  Fitzgerald's  [ 4 ] definition  of  comparabi 1 ity  of  damage 

defines  a natural  order  on  £,  ,,  the  set  of  relaxed  states  in  the 

1 rel 

<\y 

unity  configuration. 

Results  of  experiments  by  Penny  [ 6 1 motivate  the  discussion 
of  residual  stress  and  permenent  set.  Finally,  it  is  shown  how 
Noll's  theory  of  semi -elastic  material  elements  can  be  readily 
extended  to  permanent  memory  materials  with  semi-elastic  range,  PMSE. 


3.2  Bodies,  Processes  and  States 


Let  ^ be  any  set,  A function  P of  the  type 
P:[0,dp]-@.  dpcR+ 

will  be  called  a process  with  values  in  • Tbe  number  dp  will  be 
called  the  duration  of  the  process  and  the  values 


P1  = P(0),  Pf  = P(dp) 


will  be  referred  to  as  the  initial  and  final  values  of  the  process. ( 
If  Ge^  and  teR+,  we  define  a process  G,^:  [0,t]  is 

defined  by 

G(t)(r)  = G for  re[0,t] 


and  call  it  the  freeze  of  duration  t at  G. 

Let  P-j  and  P2  be  two  processes  such  that  p|  = P2-  A new  Pr°cess 

P.*P„  of  duration  d„  + dn  is  then  defined  by 
C-  Pi  Po 


(P1*p2)(t) 


P-| (t)  if  te[0,dp  ] 

P2(t-dp))  if  te[dpi.dpi  + 4pl. 


(3-1) 


p 1 *P 2 is  called  the  continuation  of  P-j  with  P2-  Continuation  is 
an  associative  process,  i.e.. 
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0VP2)*P3  = P-,*(P2-P3) 
if 

p(  = p2  = (P2*P3)'  and  - (P1*P2)f  = P3. 


Processes  of  duration  zero  are  not  excluded.  Of 
course,  all  such  processes  are  of  the  form  G^,  i,e.,  they  are 
freezes  of  duration  zero.  It  is  clear  that  for  any  process  P, 


p ■ P(0)*P  ■ P*P(0) 


Definition  1 : A body  element  is  a triple  , H ) where  0" 

is  a finite-dimensional  real  vector  space,  ^ is  a closed  and  connected 
subset  of  Sym+(,^^/**;  and  n is  a class  of  processes  with  values 
in^  which  satisfies  the  following  conditions: 


(PI)  Any  freeze  at  any  Ge^  belongs  to  H . 

(P2)  If  P belongs  to  IT  , so  does  every  segment  of  P. 

(P3)  n is  closed  under  continuation,  i.e.,  if  P-j^^H  and 
p[  = Pj,  then  P1*P2en  . 

(P4)  Any  two  elements  Gj  ,G2 can  be  connected  by  a process 
in  II.,  i.e,,  there  is  at  least  one  Pen  such  that  G-j  = P1  and 


When  speaking  about  a body  element  sr.  it  is  understood 
that  f/  is  endowed  with  a structure  defined  by  the  prescription  of 
^ and  n . The  elements  of  will  be  called  the  configuration  of 
Pj  and  the  processes  in  II  the  deformation  processes  for  Qj  . 
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If  it  is  assumed  that  $ is  arcwise  connected,  one  such  classn  that 
satisfies  (PI)  - (P4)  is  the  class  of  all  processes  P that  are 
continuous  and  have  a piecewise  continuous  derivative  P in  the  sense 
that  P(t)  exist  for  all  but  a finite  number  of  te[0,dp]  and  P has 
left  and  right  limits  for  all  te  [0,dp]. 

Most  of  the  development  herein  utilizes  the  concept  of  intrinsic 
stress.  This  describes  the  contact  forces  acting  on  an  element 
intrinsically,  that  is,  without  any  reference  to  a frame  space. 
Objectivity  is  thus  automatically  satisfied  throughout.  Full  details 
may  be  had  in  Noll  [ 5].  However,  what  has  been  said  here  will 
suffice  for  the  reader  to  understand  what  follows. 

The  primitive  notion  of  state  space  will  be  used  to  fully 
describe  the  condition,  configuration,  and  response  to  deformation 
processes  of  a material  body  element.  Permanent  memory  effects,  i.e., 
damage,  will  also  be  fully  described  by  the  state  of  the  material 
element. 

It  is  now  necessary  to  generalize  Moll's  definition  of  material 
element  to  a permanent  memory  (PM)  material  element. 

Definition  2:  A PM  material  element  is  an  octuple  («^S§>n,Z,G,S,p,D) 

which  satisfies  Axioms  I to  VI  below.  The  objects  of  the  octuple  are: 

a)  (^T^.n)  is  the  underlying  body  element.  Definition  1,  of 

the  material  element. 

b)  £ is  a set,  called  the  stats  space  of  the  material 

element. 

c)  G is  a mapping  G:  £ from  the  state  space  £ into  the 

configuration  space ^ . 
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d)  S is  a mapping  S:  £ -►  from  the  state  space  £ into 

the  stress  space  t?  - Sym  \^J  - ) • 

e)  p is  a mapping  p:  (£  x U)^  ->  £ , whose  domain 


(£  x = {(o,P)|ae  £,  PelT,  = G(a)} 


(3-2) 


J 


is  the  set  of  all  process  pairs  such  that  the  state  "fits" 
the  initial  configuration  of  the  process, 
f)  D is  a mapping  D:  2 from  the  state  space  £ onto 

the  damage  space  Q)  C Sym  )• 

The  elements  o of  £ are  the  possible  states  of  sr.  The  val 


ue 


G(cr)  is  the  configuration  and  S (a)  is  the  (intrinsic)  stress 
determined  by  the  state  cr.  The  function  p is  the  evolution  function 
of  ^Twhose  value  p(a,P)  is  the  state  reached  by  the  element  if, 
starting  from  the  state  cr,  it  is  subjected  to  the  deformation 

A 

process  P.  The  value  D(a)  is  the  damage  associated  with  the  state  a. 
In  addition,  let  us  now  define  after  Noll  [ 5 ]: 

• = (PeRIp1  = G}  as  the  set  of  all  deformation  processes 
beginning  at  the  configuration  Ge^. 

A 

• £g  = (ae£|G(cr)  = G } as  the  set  of  all  states  that  fit 
the  configuration  Ge^.  £q  is  called  the  G section  of  £. 

• The  response  functional  ^ of  the  element  is  defined  by 


£=SoS:  (lxH)fit^. 


If  oe£  and  then  £l(a,P)  1S  the  stress  produced  by 

the  process  P when 'the  initial  state  was  cr. 
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Ilg  and  IG  are  disjoint  for  different  values  of  G, 


moreover. 


H = U Ur 


GeS 


1 = ^ 


Also,  (E  x n)^t  is  the  disjoint  union  of  x ITg,  Ge 

The  foregoing  definitions  establish  the  state  a of  the  pm  material 
element  as  a primitive  concept,  embodying  within  it  the  configura- 
tion, intrinsic  stress,  and  damage  of  that  state.  It  is  a fundamental 

property  of  mappings,  that  a point  in  the  domain  is  carried  to  only 

/v  A A 

one  point  in  the  range.  That  is,  if  A is  a mapping,  and  a = x^A(a)  = A(t) 
or  equivalently:  with  a,  t elements  of  the  domain  of 
A(o)  f A(t)  <=»o  t x. 

Thus,  applying  this  observation  to  the  mappings  presented  in 
Definition  2,  we  find  that  we  cannot  have  differing  stresses,  configura- 
tions, or  damages  associated  with  the  same  stress  state. 

The  first  two  axioms  of  Noll  [5]  flow  directly  from  the 
intuitive  meaning  of  the  evolution  function  p.  The  first  states 
that  the  state  reached  after  a process  must  fit  the  final  configura- 
tion. The  second  expresses  the  fact  that  when  the  element  is 
subjected  to  a process  P1  and  then  to  a process  P2,  it  must  reach 
the  same  state  as  when  it  is  subjected  to  the  continuation  P-j*P2  of 

P1  P2 ' 

Axiom  I:  For  all  (a,P)e(E  x 


p(o,P)e  Epf , i.e,,  G(p(a,P))  ~ P ( 3-3 ) 


' t NxjUHfrS-..-1 


It  follows  from  Axiom  I that  G:  E ^ must  be  surjective 
because  by  condition  (P4)  of  Definition  1,  for  any  azZ 
and  any  there  must  be  process  PcK  such  that 

(cr,P)e(E  x n)fn  and  Pf  = G. 

Axiom  II:  If  P-j,P2  > creEpi  and  p|  = P2,  then 

p(o,(P-]*P2)  = p(p(a,P-j ) ,P2)  • ^3'4^ 

The  immediate  consequence  of  this  from  the  definition  of 
the  response  functional  above,  is 

£}  (a,P-)*P2)  = $ (p(a,P-j  ),P2)  • (3-5) 


The  next  axiom  is  as  basic  to  PM  materials  as  it  is  to 
Noll's  semi-elastic  materials.  It  expresses  the  assumption  that 
there  must  be  some  operational  way  to  distinguish  between  states. 
That  is,  if  two  states  are  different  but  fit  the  same  configuration, 
there  must  be  some  process  starting  at  these  states,  which  produces 
different  stress  responses. 

Axiom  III:  For  any  configuration  Ge^  , if  the  states 

o-j  ,o2  eEg  and  the  stress  response 

g(arP)  = £j(°2’P)  (3-6) 


for  every  deformation  process  PcH^,  then  the  states  g-j  = g2- 
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The  topology  for  PM  material  elements  is  now  introduced  following 
Noll  [ ].  The  natural  uniformity  of  the  stress  space  is  utilized 

to  develop  this  topology. 

Definition  3:  LetfJ^be  a PM  material  element  with  response  functional 

Z . For  every we  call  the  natural  uniformity  of 
of  the  <Sii  section,  the  coarsest  uniformity  on  which 
renders  the  mappings 


uniformly  continuous  for  all'PtT^.  . We  call  the 

natural  topology  of  the  state  space  Z - U ^2^  \ 

the  sum  of  the  topologies  on  the  sections  Z induced  by 


the  natural  uniformities  of  these  sections. 

The  natural  topology  onZ  is  such  that  all  sections^ are 
both  open  and  closed  inZ  • Hence  * Z-*  ^ is  trivially  continuous, 
because  it  is  constant  on  each  section: 

The  following  result  from  Noll  [ ] characterises  the  natural  topo- 

logy of  . 

r-  T- 

Proposition  1:  A net  ini-  converges  tO'S'^Z  if  and  only  if  S'-  belongs 

_____________________  2> 

eventually  to  some  fixed  section 

+ V Pc-TT& 

The  following  postulate  is  made  to  avoid  mathematical  complications. 
Axiom  IV:  The  sections!^, G , with  the  natural  uniformity  discribed 

above,  are  complete  spaces. 

If  a material  body  element  is  in  some  state  a and  is  then  frozen 
in  the  corresponding  configuration  G(a),  eventually  a "relaxed"  state 
will  be  reached. 

The  next  Axiom  of  Noll's  expresses  the  existence  of  this  state. 


Axiom  V:  For  all  ad,  the  limit 


lim  p(a  G(c)^j)  - *(o) 

t->33 


(3-7) 


exists  (and  is  unique  by  Proposition  2). 

The  mapping  x : E Z determinedly  Axiom  V will  be  called  the 

state  relaxation  mapping  with  the  members  of  its  range 


=Tel  ■ *<£) 


being  called  relaxed  states. 

Clearly,  G (X(o))  - 6(o),  and  since  G is  surjective,  the 

restriction  Glr  of  6 to  is  also  surjective. 

| re  1 

. A A 

If  Proposition  1 is  applied  to  the*  net  t ■>  p(o,G(o)^tj)  and 
equation  (5)  observed,  it  is  seen  that  (7)  is  equivalent  to 

lim  $(o,G(t)*P)  = £(Mo),P) 

t-X»  ' ' 

A 

for  all  Pen-,  where  G = G(a), 

To  simplify  notation  in  the  following  sections  the  following 
definitions  are  required: 

Definition  4:  The  set  of  all  relaxed  states  with  damage,  D, 

will  be  denoted  by  E , Q 

Therefore 
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EG  rel  D “ ^aeI  1 oeIrel  and  ^(o)  = D ] ■ 

Definition  5:  The  set  of  all  relaxed  states  in  the  configuration 

G will  be  denoted  by  , . 

G rel 

Therefore, 

ZG  rel  = £oeE  |g(o)  = G and  GeErelV 

Thus,  the  quantity  2^  ^ in  the  following  section  is  the  set 
of  relaxed  states  in  the  unity  configuration 


> 


3 
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3.3  Partial  Order  on  the  Damage  Space 
The  concept  of  permanent  memory  or  damage  is  embodied  in  the 
idea  that  permanent  material  property  changes  affect  the  state  of 
a material.  These  changes  are  assumed  to  occur  as  a result  of  the 
past  history  of  deformation  of  the  material  element.  Primary 
causes  of  damage  are  the  magnitude  of  the  maximum  deformation  gradient 
[ ] ] the  duration  of  the  deformation  process  at  the  maximum,  the 
temperature  at  which  the  maximum  occurs  [sl»  and,  if  applicable, 
a slow  (relative  to  the  deformation  processes)  healing  which  can 
decrease  or,  in  the  limit,  eliminate  the  damage  effects  [2,  3 ]. 
However,  we  are  not  here  limited  to  effects  which  weaken  the 


material.  Permanent  material  property  change  can  also  be  effected 
by  processes  which  increase  the  material  modulus.  Processes  embraced 
by  this  study  thus  include  aging,  cold  rolling  and  fiber  orientation. 

In  short,  the  concept  of  damage  embodies  all  permanent  material 
changes  that  are  brought  about  by  isothermal  strain  histories. 

The  physical  idea  of  damage  will  be  incorporated  into  the  math- 
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The  treatment  is  essentially  that  of  Fitzgerald  [ 4].  The  partial  order 
on  oD  will  be  defined  in  terms  of  the  natural  partial  order  on  the  stress 
space  of  . 

The  stress  space  can  be  regarded  as  the  set  of  self-adjoint 
bounded  linear  operations  on  the  three  dimensional  real  vector  space. 

We  can  thus  apply  the  treatment  of  Brown  and  Page  [ 9 ] to  develop  its 
partial  order. 

Definition  6:  An  operator  Te  $ is  said  to  be  positive  if  and 

only  if  (Tn,  n)  >0 

for  all  vectors  n. 

We  shall  denote  by  c^+  the  set  of  all  positive  operators  in  . 
Definition  7:  Given  T,  S e $ we  shall  write  T <_  S 

if  and  only  if  S - T e £S+- 
Obviously  T < S if  and  only  if  (Tn,  n)  < (Sn,  n) 
for  all  vectors  n [ 9 ] . 

The  stress  space  $ , is  thus  partially  ordered,  and  we  shall  use 
this  to  define  a partial  order  on  Eg  rgl . 

Definition  8:  Given  o,  t e Eg  -|  we  say  a<  t if  and  only 

if  >>^(t,P)  for  P £ "g- 

Proposition  2:  The  set  Eg  rgl  is  partially  ordered. 

Proof:  Reflexivity:  Choose  o e Eg  for  a particular  process 

we  have  ^(a,  P)  as  the  stress  response. 

By  the  reflexive  property  on 

^3  < p) 

The  process  Pe  t^G  is  arbitrary,  therefore  o < o by  definition. 
Antisymmetry:  Choose  a,  te  Eg  re^  such  that  o < t and  t < o 
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Therefore  (or  , ^ ( T , £ ^ ^(/^,  ^ ^ tS  (&>  £) 

both  equations  holding  for  all  V €.  TT  ^ . 

Thus  for  a particular  r t TT& 

SfeiV^r.E) 

by  the  antisymmetry  of  . 

Since  f1  e TT^  is  arbitrary,  we  thus  have  a = <Tby  Axiom  III. 
Transitivity:*  Choose  a.x.y  e Eg  re1  with  a < t and  r < y. 

Therefore  (.if,  5^)  CkV.  Jk  akr,2).<S(H,2^) 

both  equations  holding  for  all  Thus  for  a particular  process 

P we  have  by  the  transitivity  of  , 

SO.eV^C  4,3^ 

Since  PfcT^is  arbitrary  we  thus  have  by  definition^  N 

or  6 Tr 

Thus  ^ v<,\  is  partially  ordered. 

We  now  enlist  the  aid  of  the  following  axiom,  which  is  a generalisa- 
tion of  Nolls  definition  of  Semi  - Elastic  elements. 

Axiom  VI:  The  mapping  xVj-  is  bijective. 

Becasue  of  the  partial  order  on^vey  , this  axiom  induces  a 
natural  partial  order  on  c& 

Definition  9:  Given  o',  'Zr  £.  v-ai 

^ ancI  only  if  <s"^  Tr  . 

Thus  5*  is  partially  ordered  by  each  bijective  mapping  ^ 

The  only  step  left  is  to  make  each  partial  order  on  -^3  compatible. 

This  is  accomplished  by  placing  a restriction  on  the  choice  of  the 
mapping  D. 
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Restriction:  If  Df  ,0^  are  such  that  there  exists  ,cr,^ ^^and 

«'»  r'L- -3  with  - $ l‘A  - S ct\  Ti  Cs ' \ --  ki'  \ 

Then  if  ^ ^ C^^ 

we  must  have  C'5 ’ ) ^ i t'  \ 

Proposition  3:  The  damage  spacers  partially  ordered. 

Proof:  Reflexivity:  Choose  Then,  for  an  arbitrary  choice  of 

G*  ^ , there  exists  5 eZ  . such  that  = . 

<=>  v*iX 

Now,«" - ‘3"  . 

Therefore  TD(i~^  - (_ ^\  and  ^ 

Antisymmetry:  Choose  D-j  ,D^  € & , such  that  and  *-  . For 

an  arbitrary  choice  ofGfc^  , there  exist  states”,/?! t ZCv(t^uch 

A A 

that  D('-m)  = and  D(^)  = D^-  Therefore  by  hypothesis, 
kvS^  i CXv-  (A.  ^ ^ S C^x  ^ 

Thus  and  ^ , which  implies  that  <s\  “ ^x 

therefore  T\- 

Transitivity:  Choose  , such  that  ^ and  ^ D^. 

For  an  arbitrary  choice  ofGit  there  exist  states^  ,<3^,^^ 

, such  that  D^i)  = , D(s/)  = and  D(s\,)  = D^. 

We  have  that 

3>Ot\  * S C<3^\  cx^cX 
Therefore 

<=r>  r^x  ~^a> 

giving  that  S' , Thus  D-j  £ D^. 

We  can  see  that  a natural  consequence  of  the  above  is  that  if  D-j  £ D?, 
••on  there  exists  a configuration  G and  states  or,  e ve\  , such  that, 

^ A 

D|  and  D^l)  = D,,.  Then,  since v>  , , we  have  that 

3 <3?  vO  * V i?  fe  ^fc. 
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A subset  of  a partially  ordered  set-.&  is  said  to  be  a chain 
if  for  any  two  elements  D-j , one  has  either  D-j  < 0^  or  < D1 . 

Thus  the  restriction  of  the  partial  ordering  relation  to  yields 
a total  ordering. 

Elements  of  a chain  in  the  damage  spaced  are  said  to  be 
comparable.  Two  damage  states  not  elements  of  the  same  chain,  are 
said  to  be  incomparible. 

Definition  10.:  For  anyOT  ^ the  degree  of  damage 

ttrN,  if  and  only  if  ^ (J3~  ^ 

and 

for  at  least  one  ^fcV^and  one  direction  n. 

Lemma : If  «,J't  tIG  , then  D(«)  = D(2$  if  and  only  if 

^0.2}  = $0,^  V V>  e TV& 

Proof:  If  ”[)(«)  = D(2)  then  by  the  reflexivity  of  the  partial  order 

onJJ  we  have  D(>^)  £ D{Z)  and  D{Z)  f DfT) 
hence  ^ ^ 

and  ^ iSO.  £i>  V ^&TVt, 

For  a particular  process  P,  both  equations  hold,  therefore  by  the 
antisymmetric  property  of  the  partial  order  on 

V VL  tTVs 

The  converse  follows  by  symmetry. 

We  may  note  here,  in  order  to  emphasize  the  compatabi 1 i ty  of  this 
treatment  with  that  of  Fitzgerald [4]  , that 

*50,2^  V SclVs, 

is  equivalent  to  his 

< 'Sfe, ■£>,->  r <’SC't)$>W>  TJ  V tTT„ 
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and  for  all  vectors  n. 

The  proof  is  obvious  from  a restatement  of  the  latter  equation: 

^ C-j  ^ C ~ i ^ V)  ^ ^ ^ " o. 

Lemma:  The  setT^^  has  at  most  one  element  if  and  only  if 

the  mapping'  1S  injective. 

Proof:  Let  with  "g(^)  = G(£)  = G 

V e.w  l.  e <=»  j ''t  Cb  T ^ v ^ ^ 
o <3T  = T by  hypothesis. 

Let  ^ r6^G  rel  D. 

A. 

Since,  G(^)  = G (?)  it  follows,  by  the  injectivity  of 
\lvex  that  <3-=  T. 

a r 

Lemma:  If  A and  B are  relations  (not  necessarily  functions)  with  the 

same  domain  X,  thenAjy^is  injective  if  and  only  if 

A 

5>\-*^is  injective,  where 

y = \ > t X \ A\ 

- \ X l = 

Proof:  Let  O',  T fe  Xft  with  B(<r)  = 1i(2)  = B. 

hence  > ZT  t Xg  and  A^)  = A(<T)  = A 

hence  <$  = 2*  by  hypothesis. 

hence  %VAis  injective.  This  is  a reciprocity  lemma. 

The  converse  follows  by  symmetry. 

Proposition  4:  The  set of  relaxed  states  with  a given  con- 

figuration G and  damage  D,  has  at  most  one  element, 
which  is  unique. 

Proof:  The  proof  of  the  first  part  follows  from  both  lemmas 

and  Axiom  VI.  Uniqueness  follows  from  the  observation 
that,  if  either 

/G(«')  t\t)  or/D(<?)  f o(tf,  then‘s  f T. 
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We  now  wish  to  axiomatize  the  fact  that  every  relaxed  state  in  the 
unity  configuration  must  be  comparible  to  some  initial  material  state. 
The  treatment  follows  that  of  Penny  [ 6l. 

Axiom  VII:  The  elements  <sr e form  a space  of  chains  which 

are  not  everywhere  disjoint.  In  particular,  every 
state  in  2.^  is  comparible  to  a preferred  relaxed  state 
/\D  , called  the  virgin  state  of  the  PMM  element. 

We  normalise  the  damage  associated  with  the  virgin  state  X^to 
That  is,  = X- 

Axiom  VIII:  For  each  chain  X vaythere  exist  states<-T_  ,2^  such 

that  <3T  4 j 'o'*  4 G~  V G~  e'  X 

Proposition  5:  There  exist  states  and  22  such  that 

6~  -<  <3  * ji  rr  H ^ ^ 

Proof:  The  proof  follows  from  the  application  of  Zorn's  lemma. 

The  maximal  elements*  and  the  minimal  element  72  are  comparible  to 
e^ch  other  and  also  to  every  element  in  vaX.-  The  space  2 x 
can  thus  be  pictured  as  being  a space  of  chains  originating  at 
and  terminating  at  eithercs*  or  7*. 

/A 

If  we  distinguish  between  X,  ~ DC6*)  and  D(<T)  ^ X’  by 

saying  in  the  former  case  that  the  state  is  damaged,  and  in  the 
latter  case,  that  the  state  exhibits  strain  history  hardening. 

The  following  proposition  shows  that  a damaged  relaxed  state 
in  the  unity  configuration  cannot  have  tensile  residual  stresses. 
Proposition  6:  If®"  c'S-i  then 

a)  l<D(o)=$  S(o)  < 0 

b)  D(o)  < 1 ^ 0«  S (o) 
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Proof:  To  prove  a)  we  compare  a with  XQ  e .j  Since  D(Xq)  = 1, 

we  have  D(xq)  < D (a).  Thus,  by  definition,  XQ  < , and  so 

also  S(o,P)  < S^o’  P)  for  P LnG'  Up0n  restr1ct1n9  p t0 
being  a freeze  at  the  identity  configuration,  we  find  that 
S (a)  < S(XQ) 

But  since  \ is  the  virgin  state  and  thus  has  no  residual  stress, 

A 

we  find  that  S(o)  < 0 .The  case  of  initial  stress,  SQ,  follows 
readily.  The  proof  of  bj  follows  by  symmetry. 

3.4  Subclasses  of  Permanent  Memory  Materials 
Within  the  general  designation  of  Permanent  Memory  Material,  three 
main  subclasses  can  be  discerned.  These  are  Thixotropic  solids.  Thixo- 
tropic fluids  and  Permanent  Memory  Materials  with  Semi-Elastic  range 
(PMSE) . 

Thixotropic  Sol  id:  A PMM  element  is  said  to  be  a Thixotropic  solid 

•A 

element  if,  given  a e <-=g,  G(a)  = Gq  for  some  Gq. 

We  can  see  that  given  a stress  this  material  adopts  a unique  relaxed 
configuration  independent  of  the  damage. 

Proposition  7:  For  a Thixotropic  solid  element,  if  a e Ere-j  s=g, 

then  G(a)  = 1 

A 

Proof:  Since  XQ  e s_g,  and  since  G(xq)  = 1^  the  conclusion  follows. 

A thixotropic  solid,  therefore,  cannot  have  permanent  set  nor  residual 
stress.  The  damage  sustained  by  this  material  occurs  entirely  in  the 
"viscous"  phase  and  not  at  all  in  the  "structural"  phase. 

Thixotropic  Fluid:  A PM  material  element  is  a Thixotropic  fluid  element 

if  and  only  if  for  all  o e £rel  whose  configuration 
G = G(o)  satisfied  det  (G)^  = 1 

rJ  A 

where  G is  the  associated  bilinear  form  for  G,  S(o)=0 


3.5  Permanent  Memory  Semi-  Elastic  Materials  (PMSE) 

A PM  element  is  a PMSE  element  if  and  only  if  G\  is  injective. 

,S 

Proposition  8:  For  a PMSE  element,  if  <5,  Z e t ^ythen 

c=4>  i>C~) 

Proof:  Chooses-,  "C  e 21^  fi 

By  hypothesis  G\v  „is  injective,  thus 

*-  V*A  ,1 

G(<sj  = G(2^'=^><sr  = r c=q>  D(<5)  = 'd  [Z] . 

Corollary:  For  a PMSE  element,  ifS,  2“  <e  2. raA  ^ then 

D(«)  f D(Zr)«*G(or)  f G(2}. 


This  implies  that  relaxed  states  with  different  damages  have  different 
stress  free  configurations. 

-a 

Corollary:  For  a PMSE  element,  if  v<j,  with  D(<3")  f then 

S'  must  have  permanent  set. 

Proof:  This  follows  by  letting  Z = >>0in  the  above. 

Therefore 

D(<5)  1 1 -=>  G(®)  f J,. 


Proposition  9;  For  a PMSE  element,  if  o,  Zz  , 

then 

D(o)  = D(D4=>  a = 2T4=«>S(o)  = S (Z) . 

A A 

Proof:  All  has  been  proven  except,  S(a)  = S(t)>=<'  ° = Z, 

A I 

and  this  follows  from  the  injectivity  of  S \ ^ 
(see  reciprocity  Lemma  ) 


G rel 


Lemma:  For  a PMSE  element,  if  a e .j  ^ with  D f 1 then  the 

'V 

stress  free  configuration  is  not  the  identity  configuration. 
Proof:  Let  a e E.  . D with  D f 1.  Suppose  is  the  stress  free 

A/ 

configuration  of  o.  Comparing  a with  X , the  virgin  state, 
we  see  that,  a,  XQ  e Ere] ^ s = 
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Since  G(cs)  = G(A  ) = 1 , we  have  by  proposition  8 that 

O -v/ 

D(o)  = D(Xq)  = 1 

This  contradicts  the  hypothesis  of  the  proposition.  Therefore 
^ is  not  the  stress  force  configuration  of  a . 

Proposition  10:  For  a PMSE  element,  if  a e En  ^ n with  1 < D then 

A 

< S(a)n , n > < 0 

with  strict  inequality  holding  for  at  least  one 
direction  n. 

Proof:  By  Proposition  5 

S (a)  < 0. 

However,  since  1 is  not  the  stress  free  configuration , inequality 

rv 

cannot  hold  everywhere  in  the  above  equation. 

We  can  see  from  the  above  that  for  a strain  history  hardening 
material,  there  must  be  some  tensile  residual  stress. 

We  may  also  note  from  proposition  9 that  if  a fluid  is  PMSE,  its 
damage  must  be  manifested  as  a change  in  volume. 


PMSE  solid:  A PMSE  solid  is  a PMSE  material  such  that  D[ 

injective. 
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Proposition  11:  The  following  statements  are  equivalent  and  thus  all 

hold  for  a PMSE  solid. 

i)  The  partial  order  on  S3  defines  a partial  order 


on  £ 


u"  rel  s 

ii)  The  mapping  D \ sre-|  s 1S  injective 

iii)  Relaxed  states  with  equal  damage  and  stress  have 
equal  configurations. 

iv)  The  set  £ , 0 contains  at  most  one  element. 

' rel  s D 
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Proof: 


y i ) «=»>  i i ) ^ Let  the  Partia^  order  in  question  be: 

If  G",  Z < , then 

’ rel  s 

s'  >.<  d(«)  ^ D(r). 

Consider®',  ^tXre-|s  with  D(6)  = D(Z"). 

By  the  reflexive  property  of  the  partial  order  on  & 
we  have 

1)(C)  ^ D(7)  and  D(?)  i D(ff). 

Therefore®'  ^ 2^  and  Z ®\  and  thus®  = Z by  hypothesis. 

-*•  °lrrel  s is  injective- 

[ii) *=Mii)l  Let®-,  rel 

A /\ 

Since  D(s)  = D(2"),  it  follows  by  hypothesis  that  S’  = T. 
Therefore  G(s)  = G(Z). 

[iii)  *=>  iv)]  Let®,  2g.Z  re1st>.  By  hypothesis  G(«)  = G(D, 

and  thus  by  Axiom®"  = T. 

j*iv)t=s'>  i)l  We  must  show  that  s is  partially  ordered  by  ^ 
where  this  is  defined  for®,  2"<=ZyeVS,  as 


G ~<  U{Z). 

Reflexivity:  Let  j Therefore  D(®)€^  and  thus 

1>(e)  * D(®)- 

S’  ^ GT 

Antisymmetry:  Let  ®,  Z % with  ® - 7 and  T & <^r 

D(<r)  $ D(D  and  DfZO*  D(«) 

'*>  D(«)  = D(?) 

But®",  T , and  thus  by  hypothesis®  = T. 

Transitivity  Let6",  T,  £ 6 ZteA  ^ withS^  T and 
Therefore  D(«)  ~ D[Z)  and  D(?)  4 1)(*) 

D(er)  ^ D(8) 

or  ^ X 
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Therefore  s is  partially  ordered  and  consequently  statements 
i ) -Vo  i v ) are  equivalent  and  by  definition  must  hold  for 
a PMSE  solid. 

We  now  induce  a natural  partial  order  on  the  configurations  of 
relaxed  states  with  equal  stresses. 

Definition  11:  If®,  T ZViiV  ^ then 
G(<s)  * G(2)<t=*  <s  f 7. 

A 

Proposition  12:  For  a PMSE  solid  element,  the  set  G(ry<Ljc^  ) is  partially 

ordered. 

» ^ i 

Proof:  This  follows  from  the  injectivity  of  G \— 

Corollary:  If®',  ? <s  Zy<iJL  ^ then 

G(s)  £(?)**>  D(s)  £ D(D. 

The  proof  is  immediate,  thus  for  a PMSE  solid,  the  stress  free  states 
are  ordered  by  the  damage. 

We  see  therefore  that  incompressible  damaging  PMSE  solid  fluids 
cannot  exist. 

We  see  from  Propositionll  that  for  a PMSE  solid  s = o,D 

has  at  most  one  element.  The  class  of  materials  represented  by 
having  many  elements,  is  the  class  of  fluids  that  can  sustain  damage  in 
the  viscous  component  of  their  structure.  PMSE  fluids  can  exist, 
but  any  damage  sustained  must  be  manifested  as  a change  in  volume.  The 
material  classes  of  Mullins  elastic  and  inviscld  Fluids  are  those 
classes  such  that  every  state  is  relaxed.  The  material  classes  of 
SE  solids  and  SE  fluids  encompass  together  the  Noll  SE  class  of  materials. 
These  classes  are  characterised  by  their  inability  to  sustain  damage. 

3.6  Metrics  on  the  State  Space 

We  now  turn  to  a consideration  of  suitable  metrics  on  the  state 
spaced  . A metric,<^Cv^,  on  Z can  be  chosen  if  and  only  if,  the 
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stress  functional 

£(  -,£)  l -*>  $ 

is  continuous. 

Let  us  now  consider  a restriction  on  this  metric. 

Let  az(-,*)  = ar(-.-)  + aG(*.*)  + aD('’,)* 

where 

3 (•,•)  is  the  relaxation  metric, 

is  the  configuration  metric, 
b 

aD( * , * ) is  the  damage  metric 
and  where 


a) 

ar(alsa2)  = 0 

for  all  <Xj , $2e  £rei 

b) 

aD(a-|,a2)  = 0 

A A 

whenever  D(o-|)  = D(a2) 

c) 

aG(cfi  ,c?2)  = o 

A A 

whenever  G(a-j)  = G(a2)  . 

Proposition  13:  Every  relaxed  state  ae  E^  is  an  accumulation 
point  of 

(s,aE( - , - • 

Proof:  Choose  oe  • 

There  exists  a ae  E such  that 

A A ^ A 

0 = A (a)  = lim  p(a,G(a) ) . 

t~KO  '*•' 

Construct  the  net  i -*■  such  that 

A A 

01  = p(a,G(a)^-j) 

A 

. ‘ . lim  a-  = a 

1 -+CO 
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a - is  eventually  in  every  arbitrary  neighbor- 
hood of  a. 

A 

a is  an  accumulation  point  of  E. 

Proposition  14:  The  relaxation  metric  3r(-,*)  satisfies 

O /S 

ar(o,a)  = ar(a,t) 


A A A 

where  a = A(a),  for  all  oeE  and  for  all  xe  Erei . 

A A 

Proof:  Let  ae  E and  A(cr)  = a,  then 

ar(a,a)  < ar(a,x)  + a^x.a) 

= 3r(a,x) 

Also 

A A A A 

ar(a,T)  < ar(a,a)  + ar(a,x) 

= 3r(a,a) 

ar(a,a)  = ar(a,x)  for  all  teE^ 


by  a). 


by  a). 


Q.E.D. 


From  the  above  observations  it  is  seen  that  ar(-,-)  only  measures 

differences  in  relaxation  between  states. 

Some  particular  forms  for  these  metrics  will  now  be  discussed. 
First,  a preliminary  definition  from  Noll  [5  ]. 


Definition  12:  A function  H of  the  type  H:R+  will  be 

called  a history  with  values  in  Its  value 

H ( 0 ) = Hf 

will  be  called  the  final  value  of  H.  Let  give  rise  to  the  state 
Oj  and  H2  give  rise  to  the  state  cr2  in  what  follows. 

Let  us  now  define 


a ' (h1  »h2)  = aE(ora2), 

^,r(H-1,H2)  = ar(CTlfa2), 

3'G(Hi>H2^  = aG^al’a2^  and 
a q(^1»^2^  = 3G(a-j,a2) 


These  metrics  could  assume  the  following  particular  forms  [50]: 

a,D(H1,H2)  = am(HrH2)  + ae(H1  (o),h2(o)) 

+ ap(HrH2)  , p > 2 


where  the  norm 

I I H(s ) | |e  = [trH(s)2]14 , 


ae(H1(0),H2(0))  = I |H1 (0)  - H2(0)||e  , 


am(Hi,H2)  = 


ess  sup 
se[o,°°) 


11^(5)11, 


ess  sup 
se[o,®) 


I |H?(s) | | 


and 

ap(H! »h2)  = t f hP(s) | |H1 (s)  - H2(s)||eP  ds]1/P  . 


8A 


The  quantity  h(s)  is  the  fading  memory  influence  function  of  Coleman 
and  Noll . 

This  form  would  imply  that  damage  was  dependent  on  the  differ- 
ence in  the  maximums  (am),  the  difference  in  the  instantaneous 
viscosity  ( ) , and  a weighted  average  of  the  histories  (dp). 

The  following  could  also  be  chosen 


a,G(HrH2)  = ae(Hj  (o)  ,h2(o)  ) 


and 


a‘r(HrH2)  = a2(H15H2) 
= [ 


h2(s)||H1(s)  - H2(s)  | |e2  ds]*5. 


The  above  represent  possible  forms  for  the  metrics  3'  and 

b r 

d'g.  Again  it  must  be  noted  that  the  particular  form  chosen  must  be 
such  as  to  cause  the  stress  functional  to  be  continuous. 


3.7  Accessibility 

It  is  necessary  to  borrow  again  from  Noll  [ 5 ] to  develop  the 
idea  of  accessibility. 

The  set  of  all  states  that  can  be  reached  from  a state 
is 

(p(o,P)  |PElIg(a)}  . 

Let  us  now  define 

cl{S(a,P)|PE^(a)} 
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where  the  closure  is  in  the  sense  of  the  natural  topology  defined 

above.  If  a state  xg  I , then  the  sfate  t is  accessible  from 

o 

the  state  o.  Thus  x is  accessible  from  o,  if  every  neighborhood  of 
x contains  states  that  can  actually  be  reached  from  a. 

Transitivity  of  accessibility  is  equivalent  to  the  statement 

that,  if  tg  I then  £ Cz  . This  is  Noll 's  Proposition  13.1 . 

ox  a 

The  next  proposition,  Noll  [5],  is  a direct  consequence  of 
(7)  and  the  transitivity  of  accessibility. 

Proposition  15:  For  every  oe  I,  if  xg  I , then  x (x)g  z , 

' o a 

i.e.,  if  x is  accessible  from  o,  so  is  x (x).  In  particular, 
x (oJgZ^.  We  can  apply  Noll's  Proposition  10.5  to  PM  materials. 

Proposition  16:  In  a PM  material  element  two  states  of  equal 

damage  are  mutually  accessible. 

This  proposition  tells  us  about  accessibility  within  the  sets 
lD,  however,  the  question  of  accessibility  between  states  of  differing 
damages  is  as  yet  untouched.  Fitzgerald  [ 4 ] again  supplies  the 
appropriate  axiom. 

Axiom  IX:  Given  any  two  comparable  relaxed  states  o and  x then 

with  D(o)  > D(x)  we  say  that  0 is  accessible  from  x,o  gIa,  and  x is 
not  accessible  from  a , x ^ Ia.  Incomparabe  relaxed  states  are 
mutually  inaccessible,  i.e., 


o i Ia  and  x i Za 

r X ' o 


when  a and  x are  incomparable. 
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Aging  effects  are  neither  included  nor  excluded  in  this  treatment 
since  (9)orders  damage  only  in  relaxed  configurations.  See  Fitzgerald  [ 4 ] 
for  aging  considerations  on  PMSE  materials. 

3.8  CONSTITUTIVE  EQUATIONS 

The  previously  referenced  papers  [2,  3]  provided  a first  approximation 
to  a constitutive  equation  compatible  with  the  concept  of  damage  presented 
herein. 

Hufferd  [7]  has  shown  a computational  scheme  for  obtaining  explicit 
characterizations  from  experimental  data.  His  scheme  is  now  being  tested 
under  an  Air  Force  Rocket  Propulsion  Laboratory  contract  with  United 
Technology  Corporation. 

Lee  [8]  has  adapted  the  general  constitutive  equation  to  include 
transient  temperature  effects  and  has  successfully  modeled  the  results 
of  some  simultaneous  thermal  change  - straining  data  obtained  from 
Lockheed  Propulsion  Company  on  their  Air  Force  STV  program. 

Thiokol  Chemical  Corporation  in  collaboration  with  W.  L.  Hufferd 
is  currently  using  an  extension  of  [2,  3]  to  model  chemical  aging  effects. 

For  the  relaxed  states  discussed  herein,  the  fading  memory  terms  of 
[2,  3]  vanish  and  the  stress  tensor,  S,  as  a function  of  damage,  D,  and 
a strain  measure,  C,  become 

S = U(D,C) 

Equicontinuity  of  S in  D and  C with  the  stress  free  virgin  state 
then  permit  the  above  equation  to  assume  the  form 
S = U-| (0,  C)C  + u2(d,  C)D 

For  a positive  damage,  D,  the  material  element  in  a stress-free  con- 
figuration, C*,  then  reduces  the  above  equation  to 

U2(D,C*)D  = -U](D,C*)C* 


j»-V  ..  . . ' 
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From  [2]  it  can  be  shown  that  U2  has  a unique  inverse,  thus 
D = -U"1 (D,  D*)  ^ (D,C*)C* 

For  strain  histories  in  the  semi-elastic  range  of  a material  then, 

S = U-j  (D,  C)C  -U'1  (D,C*)  U] (D,C*)  U2(D,C)C* 

For  a pseudo-linear  semi-elastic  range  approximation 
U1 (D,C)  = U1(D) 

U (D-C)  = U-j  (D) 

hence,  S = U]  (D)  C - U2\d)  U^D)  U£  (D)C* 

For  isotropic  materials,  the  fourth  order  tensor  moduli  U-j , U2  have  the 
usual  scalar  decomposition  so  that  we  then  have 
S = U-j  (D)  (C-C*)  = U-j  (0)  C 

where  C = C-C*  is  effectively  the  strain  measure  taken  with  respect  to 
the  permanent  set  configuration. 

It  is  this  strain  measure  that  is  used  in  the  figures  of  Suli joadikusumo 
[l]  where  U-|(D)  shows  the  "softening"  effect  of  damage  which  accords  with 
our  definition  of  damage  or  herein. 

3.9  CONCLUSIONS 

A definition  and  ordering  of  the  damage  tensor  has  been  given 
mathematical  form.  The  definition  is  operational  in  that  laboratory 
tests  follow  naturally  from  the  definition. 

Several  sub-classes  of  materials,  within  the  Permanent  Memory  class 
are  elucidated.  Thixotropic  Solids,  Thixotropic  Fluids,  and  PMSE  materials. 

The  PMSE  class  exhibits  permanent  set  in  the  stress  free  configura- 
tion after  damage.  PMSE  fluids  exhibit  a volume  increase  after  damage. 

The  PMSE  class  exhibits  residual  stresses  in  the  idem  configuration 
after  damage.  This  residual  stress  is  nowhere  tensile. 

For  strain  histories  which  do  not  cause  the  damage  to  increase,  PMSE 
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materials  behave  as  Noll's  SE  group. 

The  PM  thixotropic  classes  exhibit  damage  through  their  "viscous 
components".  That  is,  damage  can  only  be  observed  during  a process 
departing  from  a relaxed  state,  transient  tests. 

PMSE  materials  exhibit  damage  through  changes  in  their  "structural" 
or  "elastic"  phase,  the  damage  may  be  observed  in  the  relaxed  state. 

PMSE  materials  may  also  exhibit  damage  effects  in  their  "viscous" 
phases. 

The  classes  herein  described,  in  addition  to  the  above,  include 
Mullins  elasticity  and  reduce  to  Noll's  SE  and  (non-linear  and  linear) 
elasticity  as  well  as  linear  fluids  on  the  limit. 

It  is  hoped  that  a proper  mathematical  structure  reflecting 
observed  mechanical  phenomena  has  been  herein  developed. 
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